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ABSTRACT 


This  report  presents  the  results  of  a  research  program  to 
develop  a  procedure  for  calculating  propulsion  system 
installation  losses.  These  losses  include  inlet  and  nozzle 
internal  losses  and  external  drag  losses  for  a  wide  variety 
of  subsonic  and  supersonic  aircraft  configurations  up  to 
Mach  4.5.  The  calculation  procedure,  which  was  largely 
developed  from  existing  engineering  procedures  and  experimental 
data,  is  suitable  for  preliminary  studies  of  advanced 
aircraft  configurations.  Engineering  descriptions,  equations, 
and  flow  charts  are  provided  to  help  in  adapting  the 
calculation  procedures  to  digital  computer  routines.  Many 
of  the  calculation  procedures  have  already  been  programmed 
on  the  CDC  6600  computer.  Program  listings  and  flow  charts 
are  provided  for  the  calculation  procedures  that  have  been 
programmed.  The  work  accomplished  during  the  program  is 
contained  in  four  separate  volumes.  Volume  I  contains  an 
engineering  description  of  the  calculation  procedures. 

Volume  II  it  a  programmers  manual  containing  flow  charts, 
listings,  and  subroutine  descriptions.  Volume  III  contains 
sample  calculations  and  sample  input  data.  Volume  IV 
contains  bookkeeping  definitions  and  data  correlations. 
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SECTION  I 


INTRODUCTION 


The  Propulsion  Installation  ai  d  Table  Assembly  Program 
(hereafter  called  PITAP)  is  designed  to  correct  uninstalled 
engine  data  for  propulsion  related  losses.  The  intended  use 
of  this  program  is  for  advanced  aircraft  configuration 
analyses.  The  degree  of  sophistication  is  consistent  with 
the  type  of  data  now  available  for  preliminary  studies  of 
advanced  aircraft  concepts. 

PITAP  consists  of  three  main  parts,  related  as  shown  in 
Figure  1.  Portions  of  the  PITAP  procedure  have  already 
been  programmed  for  the  CDC  6600  computer  and  are  now 
operational  on  the  Air  Force  Computer  System  at  WPAFB .  The 
calculation  procedures  contained  in  the  existing  computer 
program  are  indicated  in  Figure  1  by  a  dashed  line.  The 
remaining  procedures  exist  as  documented  methods  suitable 
for  computer  programming.  The  existing  computer  program 
is  suitable  for  calculating  installed  propulsion  system 
performance  corrected  for  the  most  frequently  encountered 
installation  losses.  The  programming  of  the  remaining 
procedures  will  improve  the  flexibility  of  the  program  and 
provide  greater  capability  for  making  tradeoff  studies 
and  more  options  in  the  types  of  input  data  required. 

Several  options  are  already  available  for  the  use  of  the 
existing  program,  depending  on  the  desires  of  the  user  and 
the  type  and  quantity  of  input  data  available.  The  various 
options  are  shown  in  Table  I. 

The  procedure  has  the  capability  to  correct  for  installation 
losses  associated  with  most  types  of  subsonic  and  supersonic 
aircraft  designed  for  Mach  Numbers  from  0  to  4.5. 
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EXISTS  NOW  IN 
OPERATIONAL  PROGRAM 


Figure  1:  PIT AP  PROCEDURE 


T*bie  /.  CALCULA  TION  PROCEDURE  OPTIONS 


SECTION  II 


SUMMARY 


An  analytical  program  was  conducted  to  develop  and  document 
a  calculation  procedure  that  can  be  used  to  correct  uninstalled 
propulsion  system  performance  data  for  installation  effects. 

The  program  consisted  of  three  phases:  (I)  Methods  Survey 
and  Data  Collection,  (II)  Formulation  of  Correction  Methods, 
and  (III)  Sample  Calculations.  During  Phase  I,  lata  and 
methods  were  collected  and  reviewed  to  determine  their 
usefulness  in  developing  a  general  correction  procedure. 

During  Phase  II,  the  correction  methods  were  formulated  and 
documented.  Also,  during  Phase  II,  the  main  part  calculation 
procedure  was  transmitted  to  the  Air  Force  and  has  been 
made  operational  on  the  WPAFB  CDC  6600  computer.  This  program 
can  be  used  to  make  installed  propulsion  system  performance 
calculations,  using  maps  of  inlet  performance  characteristics, 
tabulated  engine  data,  built-in  nozzle/afterbody  drag  routines, 
and  the  nozzle  internal  performance  supplied  by  the  engine 
manufacturer.  The  computerization  of  the  additional  calculation 
procedures  documented  during  the  study  will  expand  the  program 
capability  and  allow  greater  flexibility  in  types  of  input 
data.  It  will  also  permit  calculations  to  be  made  ol  inlet 
and  nozzle  data  where  no  experimental  data  exist,  using 
existing  engineering  methods  based  on  theoretical  and  semi- 
empirical  data  correlations.  Phase  III  consisted  of  samole 
calculations  to  demonstrate  the  usefulness  of  the  calculation 
procedure.  Two  configurations  were  used  for  these  sample 
cases:  an  F-4J  and  a  LWF  type  configuration. 

The  calculation  procedure  has  been  uevduped  in  the  form  jf 
several  modules  to  permit  new  data  and  methods  to  be  ad  -d 
as  they  become  available. 
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SECTION  III 


PROPULSION  SYSTEM  PERFORMANCE 


3.1  GENERAL  DESCRIPTION 

Installed  propulsion  system  performance  data  represents  the 
performance  of  an  engine  installed  in  a  particular  airplane 
configuration.  Inlet  and  nozzle  internal  performance  and  drag 
are  included,  and  bleed  air  and  power  are  extracted  -from 
the  engine  for  aircraft  systems.  These  represent  average 
values  for  aircraft  "steady-state"  performance  calculations. 
Performance  data  is  calculated  for  the  entire  engine-airframe 
flight  envelope  and  engine  power  settings  from  maximum  to 
idle. 

PITAP  can  be  used  in  several  ways (Table  I) to  produce  installed 
performance  data,  depending  upon  the  form  in  which  engine 
data  is  available.  An  engine  computer  program  can  be  used 
directly  with  PITAP,  or  it  can  be  used  with  maps  generated 
by  PITAP.  PITAP  can  also  be  used  with  tabular  data  obtained 
from  the  engine  computer  program;  i.e.,  each  computer  program 
is  run  separately.  Figures  2,  3  and  4  show  the  various 
options  that  are  possible  with  PITAP  and  an  engine  computer 
deck. 

If  a  computer  deck  is  available  for  an  engine,  the  simplest 
procedure  for  obtaining  installed  performance  is  to  use  maps 
generated  by  PITAP  (Figure  2) .  The  engine  main  program  can 
be  altered  to  accept  map  inputs ,  and  the  output  of  the  engine 
program  is  modified  to  include  inlet  and  nozzle  drag  in  the 
net  thrust  calculation.  The  engine  deck  main  program  can 
be  altered  to  include  a  table  assembly  capability. 

Figure  3  shows  the  relationship  between  PITAP  and  the  engine 
program  where  either  program  can  be  used  as  a  subroutine  to 
produce  installed  data.  Figure  4  shows  PITAP  being  used  as 
an  independent  program  with  a  specified  input  format  for 
engine  data. 

Several  possibilities  for  use  of  PITAP  have  been  described; 
however,  there  are  really  only  two  basic  methods  of  using 
PITAP  that  affect  the  PITAP  structure:  (1)  Using  PITAP  with 
tabulated  engine  performance  data  as  input,  and  (2)  Using 
PITAP  with  an  engine  cycle  match  deck  providing  engine 
performance  data. 


INLET 
RECOVERY 
MAP 


INLET 

DRAG 

MAP 


Figure  3.  PI  TAP  PROCEDURE  USING  PROGRAM  SUBROUTINES 


Figure  4:  PITAP  PROCEDURE  AS  AN  INDEPENDENT  PROGRAM 


The  engine  data  required  for  the  tabulated  version  of  the 
installation  program  are  thrust,  fuel  flow,  airflow,  nozzle 
pressure  ratio,  throat  and  exit  area  and  thrust  coefficient. 

The  installation  program  is  designed  for  engine  inputs  with 
MIL  5008B  inlet  recovery  and  corrected  mass  flow.  PITAP 
calculates  inlet  recovery  and  drag,  and  corrects  the  engine 
data  for  the  desired  recovery.  The  exhaust  pressures  and 
areas  are  required  for  nozzle  drag  calculation.  The  nozzle 
thrust  coefficient  is  required  if  the  nozzle  performance 
is  to  be  changed. 

The  tabulated  version  of  the  PITAP  installation  program 
can  also  be  used  if  only  thrust,  fuel  flow  and  air  flow  are 
available.  Methods  are  included  in  the  PITAP  procedure  to 
calculate  nozzle  pressure  ratio,  throat  and  exit  areas. 

These  methods  are  described  in  Section  "*.4.  The  calculated 
exit  areas  and  pressure  ratio  will  be  in  error  because  the 
value  of  the  exhaust  nozzle  coefficient  is  not  known;  however, 
the  resultant  error  in  nozzle  drag  calculation  is  still  within 
the  accuracy  required  for  preliminary  design  performance 
calculations . 

If  PITAP  is  used  with  an  engine  computer  program,  the 
subroutines  of  PITAP  which  are  used  to  make  inlet  recovery 
change  calculations,  nozzle  calculations,  and  thermodynamic 
calculations  are  not  required.  This  results  ir.  a  smaller 
program  with  different  logic. 

Either  version  of  the  present  programs  can  be  used  for  installed 
performance  calculation  on  turbojets,  mixed  flow  and  separate 
flow  turbofans.  Scrubbing  drag  calculation  procedures  are 
available  as  a  separate  program. 

A  more  detailed  discussion  of  the  options  available  using  the 
PITAP  program  is  contained  in  Sections  3.2  and  3.4. 

3.2  CALCULATIONS  USING  TABULATED  ENGINE  DATA 

This  is  essentially  the  existing  version  of  the  PITAP 
program  that  is  operational  on  the  WPAFB  computer.  Tabular 
inlet  input  data  are  used  to  represent  a  series  of  inlet 
performance  plots.  The  complete  version  of  PITAP  will 
eventually  include  procedures  for  calculating  inlet  performance 
plots.  Generalized  methods  are  included  in  the  program  to 
calculate  C-D  nozzle  afterbody  drag.  The  nozzle  internal 
performance  is  the  same  as  provided  by  the  engine  manufacturer. 
The  required  engine  data  is  in  the  format  of  tables  of  thrust, 
fuel  flow,  corrected  airflow  nozzle  total  or  static 
pressure  ratio,  throat  area,  exit  area  or  velocity  coefficient 
for  each  altitude,  Mach  number  and  power  setting. 


The  engine  section  of  the  PITAP  program  calculates  only  the 
changes  in  internal  performance  due  to  changes  in  inlet 
recovery.  Changes  in  inlet  recovery  produce  a  directly 
proportional  change  in  nozzle  pressure  ratio,  airflow,  and 
fuel  flow  because  the  nozzle  throat  area  does  not  change. 
Furthermore;  it  is  assumed  that  engine  data  is  calculated 
with  MIL  STD  5008B  recovery  and  all  inlet  recovery  changes 
are  made  relative  to  that  value.  Thermodynamic  data  from 
Keenan  and  Kaye  tables  has  been  "curve-fitted,"  and 
subroutines  are  provided  to  calculate  the  thermodynamic 
properties  of  the  exhaust  gases. 


The  calculation  procedure  is  as  follows:  for  each  altitude, 
Mach  number  and  power  setting,  the  net  thrust  (?N) »  fuel 

flow  (WF)»  corrected  airflow  (W  )  •  nozzle  throat  area 

(Ag),  nozzle  exit  area  (Ag)  and  nozzle  thrust  coefficient  (Cy) 
are  given. 


Standard  atmosphere  and  MIL  Standard  5008B  inlet  recovery 
are  used  to  calculate  the  airflow  at  the  engine  face  and 
gross  thrust  is  calculated  for  the  given  engine  data  before 
any  changes  in  inlet  recovery. 


W0V 

‘  oo 


F  + 
N 


The  desired  inlet  recovery  is  obtained  from  the  inlet 
subprogram  and  the  engine  gross  thrust  (F  )  is  first 
calculated  with  MIL  Standard  recovery  gl  and  then  with 
the  calculated  recovery  (F  ) ,  the  new  value  of  gross  thrust 
is  then  found  by  ratio  g2 


F 

gNEW 


F 

gOLD 


The  ratio  procedure  is  used  to  minimize  any  inaccuracies 
that  may  be  caused  by  assuming  burner  efficiency  (nB)  varies 
only  with  the  fuel-air  ratio  for  all  engines. 

The  net  thrust  and  fuel  flow  after  correction  for  inlet 
recovery  are: 


* 


* 


J 


3 
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system  thrust  and  SFC: 

°INLET  "  dnoz  +  ^oz--,., 

REF 


describe  the  methods  used  to  calculate 
the  thermodynamic  properties,  exhaust  gas  properties,  nozzle 
gross  thrust,  nozzle  pressure  ratio,  and  inlet  recovery 
correction.  These  calculation  methods  are  required  only  for 
the  tabulated  engine  data  input  version  of  PITAP. 

3.2.1  Thermodynamic  Properties 

Thermodynamic  properties  required  for  throat  calculations 
are  obtained  using  the  functions  shown  in  Table  II.  The 
functions  listed  in  Table  II  and  which  appear  in  Sections 
3.3  and  3.4  are  "curve-fits"  of  Keenan  and  Kaye  data  (Reference 
1) .  The  gas  tables  are  primarily  used  to  calculate  exhaust 
nozzle  static  pressures  and  jet  velocities. 


and  the  installed  propulsion 


N 


A 


H 


R 


WF 

SFC*  =  R 
A  — 

NA 


Sections  3.2.1  through  3.2.5 


TABLE  II 


THERMODYNAMIC  SUBROUTINES 

H  =  HOFT  (T , FOA)  Enthalpy  as  a  function  of  temperature 

(degrees  R)  and  fuel-air  ratio 

T  =  TOFH  (H,FOA)  Temperature  as  a  function  of  enthalpy  and 

fuel-air  ratio 

PR  =  PROFH  (H,FOA)  Relative  pressure,  (Pr)  as  a  function  of 

enthalpy  and  fuel/air  ratio 

H  =  HOFPR  (PR, FOA)  Enthalpy  as  a  function  of  relative  pressure 

and  fuel-air  ratio 

C  =  COFH  (H,FOA)  Sonic  velocity  as  a  function  of  total  enthalpy 

and  fuel-air  ratio 


C 


COFHS  (H,FOA)  Sonic  velocity  as  a  function  of  static 
enthalpy  and  fuel/air  ratio 
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3.2.2  Energy  Balance  for  Exhaust  Gas  Calculations 

If  the  temperature  at  the  engine  compressor  face,  airflow, 
pressure  ratio  and  fuel  flow  are  known,  the  exhaust  gas 
enthalpy  (h)  and  relative  pressure  (pr)  can  be  calculated 
from  the  aiergy  balance: 

W2hT2  +  WfQnB  =  W18hTlg  +  W8hTg  +  WBXhTBX 

(for  either  mixed  or  non-mixed  flow  engines) 

For  mixed  flow  fans  or  a  turbojet: 

W8  =  W2  WBX  +  Wf 
(f/a)g  =  wf/(w2  "  wbx^ 

hT  =  (W2hT  +  WfQnB)/Wg  ^wBXhBX  is  considerod  negligible) 
8  2 


£  <V  f/a)a 


For  a  separate  flow  ducted  fan  (fan  nozzle  and  primary  nozzle) • 
hT  =  f  (Tt,  f/a) 2 ;  where  f/a  =  0 


^PrT^2=  f  (hT'  f/a)  2 '  where  f/a  =  0 
(Pr  'l  <pt  /PT  > 

^  Iyia  rT;  2  Tia  t2 


h  =  f  (P  ,f/a) , ;  where  f/a  =  0 
T1 8  rT  18 


W8  "  W2  “  W18  "  WBX  +  Wf 

hTg  =  (W2hT2  “  W18hT18  +  WfQnB)/W8  (WBXhTBX  ls  considei‘ed 

negligible) 


'K  \  =  f  (hm ,  f/a) 


(%) 


12 


3.2.3  Nozzle  Gross  Thrust  Calculation 


The  calculation  procedure  in  this  section  applies  to  both 
mixed  and  non-mixed  flow  nozzles. 

3.2. 3.1  Convergent  Nozzle 

The  velocity  at  the  throat  for  a  convergent  nozzle  is  a 
function  of  the  total  enthalpy  (assuming  the  throat  is 
choked'/  . 

Cg  =  f(hT,  f/a)8 

and  the  static  pressure  is  a  function  of  the  static  enthalpy 

h8  =  hT0  "  (<V 


Tg  =  f  (h,  f/a)  g 


Pr  =  f  (h,  f/a). 
r8  8 


PT  /P8  =(Pr  VPr 
T8  8  rT  8  f 


P  i  is  obtained  from  the  tabulated 
"^engine  input  data  as  f  (P.S.,  alt., 
Mm)  or  it  is  calculated  by  the  procedure 
described  in  Section  3.2.4. 


P8  PTq/  ^PTg^P8  ^ 

The  area  of  the  throat  is 

A  *=  (SSI) 

8  VPC  '  8 
and  the  thrust  is 

F,  =  V|  *  V<P8  -  Pamb> 

g 

3.2. 3.2  Convergent-Divergent  Nozzle  (fully  expanded) 

If  the  exhaust  flow  is  fully  expanded,  the  static  pressure 
of  the  nozzle  exit  is  equal  to  ambient,  and  the  exit  velocity 
is  a  function  of  the  total  to  static  enthalpy. 


%  =  Pr8  <WV 


( 
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h9  =  f  (Pr,  f/a)9 
Tg  =  f  (h,  f/ a) 9 


Since  hg  =  hg 


V9  *  2gJ  (hTg  -  V 


The  exit  area  is 

A,  =  W9  R8  I9/Pamb  V, 

and  the  gross  thrust  is 


Fg  *  W9  V9  CV 


3.2. 3. 3  Convergent-Divergent  Nozzle  (not  fully  expanded) 

If  the  exit  area  of  a  convergent-divergent  nozzle  is  less 
than  required  for  full  expansion,  the  exit  static  pressure 
will  be  higher  than  ambient.  The  throat  conditions  are 
known;  therefore,  a  guessed  exit  velocity  gives: 

hQ  =  h  -  V  2/2gJ 
y  i  g  y 

Tg  =  f  (h,  f/a) 9 


P  =  f  (h,  f/a) q 
r9  y 

P9  ‘  (P8> 

r8  1 
»9  .  A9V,  =  (PAV)9  -  (|^)9 


An  iteration  on  Vg  to  make  Wg  =  Wg  will  result  in  the  exit 
conditions  for  a  given  area. 
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Since  P9  =  P™b 

PTS/Pamb  1%)^ 

The  pressure  ratio  calculation  will  be  in  error,  an  amount 
relative  to  the  value  of  the  thrust  coefficient  (Cv  )  ,  because 
this  is  usually  unknown  if  pressure  ratios  and  exhaust 

areas  are  not  given. 
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3.2.5  Inlet  Recovery  Correction 


The  engine  corrected  airflow  remains  constant  for  any  change 
in  inlet  recovery,  and  at  any  given  power  setting,  the  nozzle 
exhaust  areas  and  burner  fuel/air  ratio  remain  constant.  The 
engine  performance  for  any  change  in  inlet  recovery  is  calculated 
by  the  following  relations: 


RATIO 

RF/RFMIL ;  where  RF/RFMIL  = 

(pt„/pt  5  /  (PT,/PT  *mil 

2  "  2  “  5008B 

W8RF 

= 

W8 

(RATIO) 

WFRF 

= 

WF 

(RATIO) 

W2RF 

= 

W2 

(RATIO) 

P8P0RF 

=3 

P8/P0  (RATIO) 

After  the  above  quantities  are  computed,  the  corrected 
quantities  (W8RF ,  WFRF ,  W2RF,  and  P8P0RF)  are  used  in  the 
subroutines  described  in  Sections  3.2.1,  3.2.2,  and  3.2.4 
to  compute  a  new  gross  thrust,  FQ  .  This  new  gross  thrust 

and  the  gross  thrust,  F_  ,  calculated  using  the  same  sub- 

G1 

routines  and  the  uncorrected  (MIL  5008B)  quantities  (W8,  WF, 
W2,  P8P0)  are  used  to  compute  a  ratio,  F_  /F„  .  This  ratio 

G2  G1 

is  then  used  as  described  in  section  3.2  to  obtain  the  new 
value  of  gross  thrust,  Fr 


3.3 


NEW 

CALCULATIONS  USING  ENGINE  CYCLE  MATCH  OECK 


This  version  of  the  PITAP  program  differs  from  the  existing 
PITAP  program  in  two  specific  areas.  The  first  is  in  the 
input/output  of  the  program  and  the  second  is  that  the 
thermodynamic  functions,  performance  and  thrust  subroutines 
are  not  required.  For  this  version  of  the  PITAP  program  using 
the  engine  cycle  match  deck,  the  PITAP  program  is  the  main 
program  and  the  engine  cycle  match  deck  is  used  as  a  subroutine. 


The  input/output  of  the  engine  program  and  PITAP  are  linked 
through  the  labeled  common  blocks  of  each  program.  The  input 


for  the 

engine  program 

requires  the 

following  data: 

PS 

= 

power 

setting  or  code 

T2 

degrees 

R  = 

inlet 

total  temperature 

P2 

lb/in*- 

= 

inlet 

total  pressure 
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Tamb 

degrees  R  = 

free-stream  static  temperature 

4 

V. 

Pamb 

lb/in2 

free-stream  static 

ETAR 

inlet  recovery 

WBFAN 

lb/sec  = 

fan  bleed  airflow 

l 

lb/sec  = 

bleed  air  flow  (high  pressure) 
extraction 

V* 

lb/sec  = 

bleed  airflow  (intermediate) 

V 

HP 

HP 

horsepower  extraction 

Power  setting,  bleed  and  horsepower  extraction  are  input,  but 
inlet  and  ambient  pressure  conditions  are  calculated  from  the 
altitude  and  Mach  input.  The  existing  version  of  the  PITAP 
Program  obtains  atmospheric  properties  from  a  subroutine  which 
provides  ambient  temperature  and  pressure  for  ISA  standard 
day  as  a  function  of  pressure  altitude.  For  non-standard  day 
conditions,  a  different  subroutine  will  be  required.  A  special 
subroutine  is  provided  with  the  engine  program  to  interface 
with  the  inlet  subroutine  to  obtain  ram  recovery. 

It  is  recommended  that  power  setting  or  power  code  inputs 
conform  with  ARP  681.B  (PS  =  100  is  maximum  augmented  thrust, 

50  is  intermediate,  20  is  idle,  etc.) 

Five  tables  are  added  to  the  current  version  of  PITAP  to 
input  fan  bleed,  LP,  IP  and  HP  compressor  bleed,  and  horsepower 
extraction  as  a  function  of  altitude  and  flight  Mach  number. 

The  following  input  variables  are  required  for  variable  cycle 
decks  or  engine  description: 


Input  Variables 


Symbol 

Unit 

Definition 

BRPDS 

- 

Design  bypass  ratio 

P3P2DS 

- 

Design  overall  compressor  pressure  ratio 

T4MAX 

degrees  F 

Maximum  turbine  inlet  temperature.. 

W/DES 

lb/sec 

Design  total  airflow 

17 


Input  Variables  (Continued) 


Symbol 

Unit 

Definition 

PCW2C 

- 

Per  cent  corrected  airflow  desired  at 
maximum  height  and  Mach  number 

ALT  MAX 

ft 

Altitude  at  maximum  flight  Mach  number 

MMach 

Maximum  Mach  number 

DASH  M 

Dash  Mach  Number 

D  ALT 

ft 

Dash  altitude 

LENG 

in 

Overall  length 

LNOZ 

in 

Nozzle  length 

DIA  INL 

in 

Inlet  flow  path  diameter 

DMAX 

in 

Maximum  engine  diameter 

WTENG 

lb 

Engine  weight 

WTNOZ 

lb 

Nozzle  weight 

The  engine  performance  data  is  run  for  a  matrix  of  Mach  numbers, 
altitude  and  power  settings. 

The  output  format  shown  in  Figure  5  uses  67  spaces  on  a 
printout  page  so  that  it  can  be  cut  to  8  1/2  x  11  size  for 
storage,  or  it  can  be  in  CRT  displays  and  remote  teletype 
systems . 

The  output  format  has  classification  headings  at  the  top  and 
bottom  of  the  page  and  is  divided  into  five  basic  groups  of 
data: 

a.  Engine  company,  engine,  and  run  identification 

b.  Engine  design  description  which  is  either  input  or  calculated 

c.  Engine  data 

d.  Installed  performance  -  inlet  and  nozzle  drag 

e.  Special  output:  i.e.,  FN/6,  SFC/  VIT  etc. 
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Figure  5:  0U  TPU  T  FORMA  t  for  PI  TAP  (C  ycl  e  ma  tch  deck  VERSION ) 


Output  Parameters 


Symbol 

Units 

Definition 

Classification  - 

engine  program  security  classification 

Engine 

Co.  Name  - 

engine  co.  name,  engine  identification 
and  run  title 

BPRDS 

- 

design  bypass  ratio 

P3P2DS 

- 

design  overall  compressor  pressure  ratio 

T4MAX 

degrees  F 

maximum  turbine  inlet  temperature 

WOOES 

lb/sec 

design  total  airflow  -  sea  level,  standard 
day 

PCW2C 

- 

percent  of  design  corrected  flow  required 
at  maximum  0^  condition 

ALTMAX 

ft 

altitude  of  maximum  Mach  number 

MMAX 

- 

maximum  flight  Mach  number 

DASHM 

- 

dash  Mach  number 

DALT 

ft 

dash  attitude 

LENG 

in 

engine  length  (including  nozzle) 

LNOZ 

in 

nozzle  length 

D1AINL 

in 

inlet  flow  path  diameter 

DMAX 

in 

maximum  engine  diameter 

WTENG 

lb 

engine  weight  (including  nozzle) 

WTNOZ 

lb 

nozzle  weight 

ALT 

ft 

altitude 

MACH 

- 

Mach  number 

DELT 

degrees  F 

temperature  deviation  from  standard  day 

Pb 

- 

power  setting 

PCODE 

- 

identification  of  power  setting  i.e., 
intermediate,  max  augmented,  etc. 
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Output  Parameters  (Continued) 


Symbol 

Units 

Definition 

FN 

lb 

net  thrust 

SFC 

lb/hr/lb 

specific  fuel  consumption 

FRAM 

lb 

ram  drag 

WFT 

lb/hr 

total  fuel  flow 

FG 

lb 

gross  tl  ~ust 

PAMB 

lb/in2 

ambient  pressure. 

P2 

lb/in2 

inlet  total  pressure 

P12 

lb/in2 

fan  discharge  total  pressure 

P22 

lb/in2 

LP  compressor  discharge  total 
(port  losses  included) 

pressure 

P2x 

lb/in2 

IP  compressor  discharge  total 
(port  losses  included) 

pressure 

P3 

lb/in2 

HP  compressor  discharge  total 
(port  losses  included) 

pressure 

TAMB 

degrees  P 

Ambient  temperature 

T2 

degrees  R 

engine  inj/et  total  temperature 

T12 

degrees  R 

fan  discharge  total  temperature 

T22 

degrees  R 

LP  compressor  discharge  total 

temper  ature 

T2X 

degrees  R 

IP  "  " 

It 

T3 

degrees  R 

HP 

M 

WBFAN 

lb/sec 

fan  bleed  air  flow 

WBLP 

lb/sec 

LP  compressor  bleed  air  flow 

WBI 

lb/sec 

IP  compressor  bleed  air  flow 

WBA 

lb/sec 

HP  compressor  bleed  air  flow 

HPX 

hp 

horsepower  extraction 

DEX 

in 

nozzle  exit  diameter 

Output  Parameters  (Continued) 


Symbols 

Units 

Definitions 

THETA 

degrees 

nozzle  boattail  angle 

ETAR 

- 

inlet  total  pressure  lecovery  ratio 

W2C 

lb/sec 

engine  corrected  air  flow 

P8P0 

- 

exhaust  nozzle  pressure  ratio 

T8 

degrees  R 

exhaust  nozzle  total  temperature 

W3 

lb/sec 

HP  compressor  air  flow 

W8 

lb/sec 

exhaust  nozzle  mass  flow 

WF4 

lb/hr 

burner  fuel  flow 

WF8 

lb/hr 

total  fuel  flow 

FGI8 

lb 

ideal  gross  thrust 

CV8 

- 

nozzle  thrust  coefficient 

DRAG 

- 

nozzle  external  drag 

AE8 

in2 

nozzle  throat  effective  area 

A8 

.  in2 

nozzle  throat  area 

AE9 

in2 

nozzle  exit  effective  area 

A9 

in2 

nozzle  exit  area 

P180PO 

- 

fan  nozzle  pressure  ratios 

T18 

degrees  R 

fan  nozzle  exhaust  temperature 

W12 

lb/sec 

fan  duct  air  flow 

W18 

lb/sec 

fan  nozzle  mass  flow 

WF18 

lb/hr 

fan  duct  reheat  fuel  flow 

FGI18 

lb 

fan  nozzle  ideal  gross  thrust 

CV18 

- 

fan  nozzle  thrust  coefficient 

AE18 

.  2 

in 

fan  nozzle  throat  effective  area 
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Output  Parameters  (Continued) 


Symbols 

Units 

Definitions 

A18 

.  2 
m 

fan  nozzle  throat  area 

AE19 

.  2 
in 

fan  nozzle  exit  effective  area 

A19 

.  2 
m 

fan  nozzle  exit  area 

FNA 

lb 

propulsion  system  net  thrust 

SFCA 

lb/hr/lb 

propulsion  system  specific  fuel  comsumption 

DINLET 

lb 

inlet  drag 

DNOZ 

lb 

nozzle  drag 

AC 

.  2 
m 

inlet  capture  area 

AOE/AC 

- 

engine  mass  flow  ratio 

AO I /AC 

- 

inlet  capture  mass  flow  ratio 

AO/AC 

- 

inlet  mass  flow  ratio 

CPSPL 

- 

coefficient  of  inlet  spillage  drag 

CDBCD 

- 

coefficient  of  inlet  bleed  drag 

CDBYP 

- 

coefficient,  cf  inlet  bypass  drag 

CDINL 

- 

coefficient  of  inlet  total  drag 

AMAX 

.  2 
in 

maximum  area  of  the  exhaust  system 

LBT 

in 

length  of  exhaust  system 

DNOZR 

lb 

reference  nozzle  drag 

CDBT 

- 

coefficient  nozzle  boattail  drag 

CDBASE 

- 

coefficient  nozzle  base  drag 

CDINT 

- 

coefficient  nozzle  interference  drag 

BETA 

degrees 

boattail  angle 

FN/DELTA 

lb 

corrected  net  thrust 

SFC/RTHETA 

.  lb/hr/lb 

corrected  specific  fuel  consumed 

Q 

lb/ft2 

free-stream  dynamic  pressure 
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3.4  OTHER  PROGRAM  OPTIONS 


3.4.1  Inlet  and  Nozzle  Performance  Maps 

A  subroutine  has  been  programmed  that  can  call  and  control 
the  existing  subroutines  required  to  generate  inlet  performance 
"maps"  from  tabulated  inlet  input  data  (representing  inlet 
performance  "plots").  The  program  listing  for  this  subroutine 
is  included  in  Volume  II.  A  subroutine  also  exists  to  generate 
nozzle  data,  using  the  internal  program  procedures,  but  the 
output  needs  to  be  revised  to  adequately  output  the  map  data 
in  the  desired  format. 

3.4.2  Engine  Manufacturers  Nozzle  Internal  Performance 

The  existing  PITAP  program  (described  m  Volume  II)  \ises 
this  procedure. 

3.4.3  Externally-Generated  Cy  Input 


This  option  to  the  basic  program  procedure  makes  it  possible 
to  drop  out  the  engine  manufacturers  nozzle  internal  performance 
(Cy)  and  substitute  a  different  Cy,  which  can  be  either  a 
constant  value  or  a  map.  v 

This  is  accomplished  by  changes  to  the  existing  Main  Program 
routine  (TEM  333)  and  subroutine  COMPUTE.  The  changes  consist 
of  the  following: 

1.  The  main  program  is  changed  to  accommodate  reading 
in  and  storing  the  new  coefficient  or  map. 

2.  Subroutine  COMPUTE  is  changed  to  call  the  stored 
coefficients  and  insert  them  instead  of  the  engine 
nozzle  coefficient. 

The  changes  are  not  extensive  and  the  logic  for  them  is 
included  in  Volume  II. 

3.4.4  Nozzle  Internal  Performance  Subroutine 

This  involves  modifying  the  Main  Program  and  Subroutine  COMPUTE 
to  eliminate  use  of  the  engine  nozzle  loss  coefficient  and 
building  in  the  logic  to  call  and  use  a  subroutine  which  has 
the  nozzle  internal  loss  calculation  procedures  built  into  it. 
The  nozzie  internal  loss  calculation  methods  are  described  in 
Section  V  of  this  document.  Engineering  logic  diagrams  are 
provided  in  the  Appendix  to  show  how  this  should  be  done. 
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SECTION  IV 


INLET  PERFORMANCE 


The  major  installation  losses  that  must  be  accounted  for  in 
correcting  uninstalled  propulsion  system  data  for  inlet  effects 
are:  (1)  total  pressure  recovery,  (2)  boundary  layer  bleed 

drag,  (3)  bypass  drag,  and  (4)  spillage  drag.  This  section 
describes  the  methods  used  to  calculate  these  corrections.  A 
general  description  of  the  existing  computerized  methods  is 
presented  first.  The  existing  computer  program  uses  input 
"maps",  representing  inlet  recovery  and  drag  as  a  function 
of  mass  flow,  to  calculate  installation  corrections. 

Following  a  description  of  the  engineering  procedures  in  the 
computer  program,  engineering  methods  are  discussed  which  can 
be  used  to  obtain  input  data  when  no  input  "maps"  are  available. 
Several  of  the  engineering  methods  have  been  programmed  for 
the  CDC  6600  computer  and  programming  flow  charts  and  listings 
for  these  are  included  in  Volume  II. 

Figure  6  presents  a  diagram  showing  how  the  inlet  calculation 
procedures  relate  to  the  installed  propulsion  system  procedures. 

4 . 1  GENERAL  DESCRIPTION 

Corrections  to  uninstalled  propulsion  system  performance  to 
account  for  inlet  recovery,  bleed  drag,  bypass  drag,  and 
spillage  drag  can  be  accomplished  using  several  different 
types  of  input  data. 

The  basic  propulsion  system  calculation  procedure  uses  "maps" 
of  total  pressure  recovery  and  inlet  drag  as  a  function  of 
engine  corrected  airflow  divided  by  capture  area  (Item  3  in 
Figure  7) .  These  maps  directly  provide  the  recovery  and 
drag  for  use  in  installed  propulsion  system  calculations . 
Corrected  airflow  is  the  parameter  used  to  match  the  inlet 
performance  to  the  engine. 

However,  most  wind  tunnel  test  data  are  obtained  in  the  form 
of  individual  ir.let  performance  "plots"  such  as  those  shown 
in  Item  2  of  Figure  7.  A  separate  subprogram  is  used  to 
convert  these  performance  plots  into  performance  maps.  This 
subprogram,  (INLTMAP)  is  described  in  Volume  II. 

If  wind  tunnel  test  data  are  not  available  for  which  to  obtain 
inlet  performance  plots,  these  plots  are  builL  up  by  theoretical 
and  semi-empirical  calculation  procedures  (Item  1  of  Figure  7) , 
that  are  described  in  this  report. 
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Figure  6:  INLET  PERFORMANCE  PROCEDURE 


COMPONENT 

PERFORMANCE 


Figure  7:  INL  ET  PERFORMANCE  EL  EM  ENTS 
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The  inlet  performance  calculations  are  performed  by  the  inlet 
subprogram,  Figure  8. 

Inlet  performance  maps  are  input  data  to  the  inlet  subprogram. 
This  subprogram  sizes  the  inlet  capture  area  (if  it  is  not 
specified)  and  converts  the  inlet  performance  maps  into  total 
pressure  recovery  and  inlet  drags  that  are  matched  to  the 
corrected  airflow  demand  of  the  engine. 

The  connecting  link  between  the  engine  data  and  the  inlet 
subprogram  is  engine-plus-secondary  corrected  airflow.  The 
sizing  routine  permits  the  inlet  to  be  sized  for  operation 
at  a  desired  inlet  mass  flow  ratio,  consistent  with  a  recovery 
associated  with  the  design  engine  corrected  airflow  demand.  A 
specified  capture  area  size  can  also  be  input  if  desired, 
instead  of  requiring  the  program  to  calculate  the  size. 

The  inlet  input  requires  twelve  tables  of  input  data  which 
describe  the  performance  characteristics  of  the  inlet.  Engi¬ 
neering  data  obtained  from  wind  tunnel  tests  and  theoretical 
calculations  are  used  to  obtain  the  inlet  performance  charac¬ 
teristics.  The  format  for  the  inlet  data  is  shown  in  Figures  9 
and  10.  Data  taken  from  these  engineering  plots  are  punched 
on  cards  for  input  into  the  inlet  subprogram. 

The  inlet  subroutine  recognizes  three  modes  of  inlet 
operation:  low-speed,  external  compression,  and  mixed 
compression.  The  low-speed  mode  is  used  only  at  very  low  Mach 
numbers,  e.g.,  take-off  conditions,  when  only  high  engine 
power  settings  are  likely  to  be  of  interest  and  inlet  drag 
is  negligible.  The  external-compression  mode  is  used  over 
the  remaining  Mach  number  regime  for  external-compression 
inlets.  It  is  also  used  for  the  remaining  subsonic  regime 
and  supersonic  Mach  numbers  up  to  the  starting  Mach  number 
for  mixed-compression  inlets.  The  mixed-compression  mode  is 
used  at  or  above  the  starting  Mach  number  for  mixed-compression 
inlets . 

a)  External-Compression  Inlets.  The  calculation  of 
recovery  and  drag  for  an  external-compression 
inlet  is  illustrated  in  Figure  9.  The  required 
performance  maps  are  input  as  tables,  as  indicated. 
Table  1  is  used  to  represent  the  effect  of  the 
airplane  flow  field  on  the  local  Mach  number 
seen  by  the  inlet.  Table  2A  gives  the  basic 
recovery /mass-flow-ratio  characteristics  of  the 
inlet.  The  minimum  Mach  number  for  which  data 
is  input  in  Table  2A  is  taken  by  the  program  to 
be  M  ,  below  which  only  the  low-speed  mode  is 
min 

used. 
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Figured:  INLET  SUBPROGRAM 
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Table  2B 
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Figure  10:  PERFORMANCE  CA  L  CULA  T/ON  FORA  MIXED-COMPRESSION  INL  c  T 
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In  the  low-speed  mode,  recovery  is  read  directly 
out  of  Table  2B  as  a  function  of  local  Mach  number 
only,  and  inlet  drag  is  neglected.  If  the  Local 
Mach  number  exceeds  ,  the  recovery  and  mass 

flow  ratio  are  determined  using  Table  2A,  Table  7 
(which  gives  the  scheduled  bypass  flow,  if  any,  as  a 
function  of  engine  mass  flow  ratio) ,  and  engine 
corrected  airflow  demand.  An  iteration  is  performed 
to  solve  simultaneously  for  the  match-point  recovery 
and  inlet  mass  flow  ratio,  as  well  as  the  engine  mass 
flow  ratio  and  scheduled  bypass  flow.  If  the  indicated 
buzz  (Table  2D)  or  distortion  (Table  2E)  limits  are 
exceeded,  an  appropriate  warning  message  will  appear, 
but  no  fatal  error  will  result.  The  bleed  mass  flow 
associated  with  the  calculated  inlet  mass  flow  ratio 
is  determined  from  Table  6A. 


After  the  required  mass  flow  ratios  are  determined, 
spillage,  bleed,  and  bypass  drags  are  found  from 
Tables  3,  4,  and  5,  respectively.  Spillage  drag 
is  the  incremental  change  in  additive  drag  and 
pressure  drag  on  the  airplane  due  to  inlet  operations 
at  mass  flow  ratios  less  than  a  reference  mass  flow 
ratio.  The  bleed  and  bypass  drags  include  door 
drags  as  well  as  momentum  losses  of  the  airflow. 

b)  Mixed-Compression  Inlets.  The  performance  calculation 
for  a  mixed-compression  inlet  is  illustrated  in 
Figure  10.  Below  the  starting  Mach  number  Mg,  the 
low  speed  mode  and  external-compression  mode  are 
used  in  the  same  way  as  in  the  case  of  an  external- 
compression  inlet.  The  mixed-compression  mode,  used 
at  or  above  Mg  ,  is  based  on  the  assumption  that  a 
closed-loop  bypass  system  is  available  to  remove 

all  excess  air.  Thus,  except  for  the  case  ot 
excessive  engine  airflow  demand,  the  inlet  mass  flow 
ratio,  bleed  flow,  and  recovery  may  all  be  scheduled 
as  a  function  of  local  Mach  number  only;  the  bypass 
system  compensates  for  changes  in  engine  airflow 
demand. 

If  the  corrected  airflow  delivered  by  the  inlet  is 
inadequate  to  meet  the  engine  demand  at  the  scheduled 
recovery,  the  program  will  permit  the  inlet  to 
operate  at  an  excessive  supercritical  margin.  The 
recovery  will  be  lowered  sufficiently  to  match  the 
engine  corrected  airflow  demand,  and  an  appropriate 
message  will  warn  the  user  of  an  undersized  inlet. 


32 


Inlet  spillage,  bleed,  and  bypass  drag  are  found 
using  Tables  3,  4  and  5  as  in  the  external- 
compression  mode.  The  data  in  the^t  tables  for 
Mach  number  equal  to  or  greater  th.>n  M„  apply  only 
for  started  inlet  operation.  ° 

4.2  TOTAL  PRESSURE  RECOVERY 

/ 

Three  separate  procedures  are  used  to  calc-.’  i.ate  inlet  total 
pressure  recovery,  depending  on  the  flight  speed  of  the 
airplane. 

4.2.1  Low  Speed 

At  takeoff  and  low  speeds  (0  5  M  i  .40),  the  primary  total 
pressure  loss  generating  mechanisms  are  sharp  lip  losses 
and  subsonic  diffuser  losses.  The  methods  used  to  calculate 
sharp  lip  losses  accounts  for  lip  bluntness,  effect  of  takeoff 
door  area,  and  inlet  flow  velocity.  The  subsonic  diffuser 
losses  are  calculated  by  using  an  input  variation  of  duct 
loss  coefficient,  e  =  AP.p/q^,  as  a  function  of  inlet  throat  Mach 

number.  The  duct  loss  variation  is  obtained  from  a  catalog  of 
data  provided  (in  Volume  IV)  for  a  range  of  configuration 
variables. 

The  sharp  lip  losses  are  combined  with  the  subsonic  diffuser 
losses  as  shown  in  the  flow  chart  of  Figure  60  to  obtain  the 
overall  inlet  total  pressure  recovery. 

The  nomenclature  used  to  calculate  the  lip  losses  for  low 
speed  inlet  operation  is  shown  in  Figure  11. 

The  required  input  quantities  for  the  tow  speed  calculation 
are  j 


lip  bluntness  parameter,  ratio  of  *ip 
radius  to  inlet  lip  hydraulic  diaiuetex  , 
dimensionless . 

engine  corrected  airflow,  lb/see. 
inlet  capture  area,  in.2 


M 


0 


local  inlet  Mach  number,  dimensionless 


A 


T 


inlet  throat  area, 


in. 


2 
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etric 

it 


'Eff 


Actual 

Effective 


(a)  SHARP  LIP 


BLUNT  LIP 


k  Low-speeo  inter  up  losscs 


AT.O. 


Takeoff  door  throat  area, 


2 

in. 


e 


subsonic  diffuser  duct  loss  coefficient, 

Ap 

e  =  T  ,  dimensionless. 


The  basic  calculation  procedure  is  the  same,  regardless  of 
the  local  Mach  number;  however,  below  M  =  .20,  inlet  throat 
Mach  number  is  based  on  the  calculated  throat  corrected 
airflow  parameter, 


VT 

T 


wT 


GEOM 


This  parameter  is  calculated  by  the  following  equation: 


W 


VT 1 

AT 

GEOM 


(1) 

Inlet  throat 
geometric  area, 


After  the  throat  corrected  airflow  parameter, 
calculated,  throat  Mach  number  based  on 
geometric  throat  area,  M  ,  is  obtained 

from-  iGEOM 


W  VT  ,  is 
F~  A 

GEOM 


W  V"e 
3  A 


.592 


_ M 

(1  +  . 2M2) 3 


Y  =  1.40 


This  equation  can  be  solved  by  computer  or  the  input  data 
curve  can  be  programmed  for  a  table  look-up  procedure.  This 
curve  is  plotted  in  Figure  12. 


Above  Mq  -  .20,  the  throat  Mach  number  based  on  geometric  area, 

M  ,  is  obtained  by  calculating  the  throat  sonic  area  ratio, 
GEOM 

(A/A*)  .  To  accomplish  this  calculation,  the  inlet  mass 

GEOM 

flow  ratio  is  first  calculated: 


(2) 
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lisaai. 


Next,  (— »)  is  calculated  as  follows: 

rGEOM 

A*>  <A  1  (A  ’  A*  q 

GEOM 


is  then  obtained  from  the  relation: 

GEOM 


A  (1  +  .2M2)3 

A*  1.728M 


y  =  1.40 


(3) 


(4) 


This  relation  is  plotted  in  Figure  13,  for  convenient  use  in 
performing  hand  calculations,  if  desired. 

The  effect  of  adding  flow  area  by  takeoff  doors  is  accounted 
for  by  adding  the  takeoff  door  geometric  throat  area  to  the 
inlet  geometric  throat  area  and  then  calculating  the  throat 
Mach  number  the  same  way  as  before. 

After  calculating  the  throat  Mach  number  based  o.\  geometric 
area,  M  ,  the  effective  throat  Mach  number,  M  ,  is 
uEOM  XEFF 

obtained  from  plotted  data  (or  table-look-up,  in  case  of  the 
computer  program)  showing  the  variation  of  M  as  a  function  Of 

XEFF 

M  and  lip  bluntness  parameter,  r/D.  An  example  of  this 

GEOM 

plotted  data  is  shown  in  Figure  14. 

The  ratio  of  effective  throat  area  to  geometric  throat  area, 

EFF  ,  is  next  obtained  from  the  following  equation: 

AT 

GEOM 


at  (—  ) 

EFF  EFF 

AT  (—  ) 

GEOM  VA* ' GEOM 

=  f  (^T 

v  EFF, 

M,p  j 

uEOM 

(5) 

Am 

Pn 

The  simplifying  assumption 

is  made  that 

T 

EFF  = 

TT. 

A,,, 

^T 

1 

GEOM 

This  assumption  is  based  on  the  following  rationale: 

If  A^  =  A2  but  there  is  a  loss  in  total  pressure  between 

point  1  and  2,  then  the  Mach  number  at  point  2  will  be 
higher  than  that  at  point  1,  because  of  the  effective 
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Figure  14:  EFFECTIVE  INLET  THROA  T  MACH  NUMBER 


reduction  in  area  ratio,  (A/A*) since  (A*) 2  has  increased.  The 

present  method  assumes  that  the  reduction  in  (A/A*)2  is  equivalent 

to  that  obtained  by  reducing  A2  and  keeping  A*  constant  between 

points  1  and  2.  Thus,  there  is  a  direct  relationship  between 
A_  and  A.  which  is  related  to  the  3  ip  loss.  This  relationship 
tFF 

can  be  expressed  as: 

?T  ^2 
T  =  *EFF 


and  since  A2  =  A^,  and  A 2  =  AT 


‘EFF 


GEOM 


at  I  r 

Thus,  the  EFF  can  be  substituted  for  T. 

A„  P„ 


GEOM 


T. 


After  throat  effective  Mach  number,  M  ,  is  obtained,  the 

*EFF 

subsonic  diffuser  duct  loss  coefficient,  t  =  AP_,/q,,  ,  is  obtained 
from  input  data  showing  the  variation  of  e  versus  M 

1  EFF 


c  is  then  used  in  the  following  equation  to  calculate  the 
total  pressure  recovery  of  the  subsonic  diffuser: 


i 


3.5 

(1  +  .2M  2) 

EFF 


(6) 


For  convenience  in  performing  hand  calculations  this  relation 
is  shown  plotted  in  Figure  15. 

Overall  total  pressure  is  next  computed  from  the  product  of 
throat  recovery  and  diffuser  recovery: 
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Figure  15:  SUBSONt 


An  assumption  of  PT  /P_,  =  .95  is  used  for  this  purpose.  The 

2  0 

procedure  operates  by  assuming  a  value  of  recovery,  calculating 
throat  Mach  number,  calculating  sharp  lip  losses,  calculating 
subsonic  diffuser  losses,  sums  the  losses  to  obtain  overall 
total  pressure  recovery  and  compares  the  calculated  value  with 
the  assumed  value.  If  the  calculated  and  assumed  values  differ 
by  more  than  .001,  a  new  value  of  P_  /P  is  assumed  that  is 

l2  lO 

equal  to  the  average  of  the  calculated  and  assumed  recoveries 
and  the  calculation  procedure  is  repeated. 

4,2,2  Subsonic  and  Transonic 

In  this  speed  range  (.4  £  M  <  1.10)  the  inlet  total  pressure 
recovery  is  determined  almost  entirely  by  the  efficiency  of 
the  subsonic  diffuser.  Boundary  layer  bleed  flow  is  usually 
not  required.  If  excess  air  is  taken  into  the  inlet,  it  is 
discharged  through  a  bypass  system.  At  Mach  numbers  below 
1.10,  spillage  is  usually  a  better  way  to  get  rid  of  excess 
inlet  airflow  because  it  has  a  lower  drag  penalty  than  bypassing 
at  these  Mach  numbers. 


The  total  pressure  recovery  calculation  procedure  is  shown  in 
the  flow  chart  of  Figure  61. 


The  procedure  involves  an  iterative  solution  for  the  diffuser 
losses  that  matches  the  airflow  through  the  inlet.  Initially, 
a  recovery  value  of  P  /p  =  .95  is  assumed  to  begin  the 

2  0 

calculation.  Using  the  assumed  recovery,  inlet  mass  flow  ratio 
(assuming  no  bypass  or  internal  bleed  is  used)  is  calculated 
by  the  following  equation: 


(w  VT\ 

rw 


<v 


.  343 


Next, 

area, 


the  throat  area  ratio  based  on  geometric  throat  flow 
,A  ,  is  calculated  from  the  following  equation: 

'a7' 


GEOM 


(Aq/Ac)  (p)  0 
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From  this  area  ratio,  the  throat  Mach  number,  M  ,  is 
obtained  from  the  relation:  XGEOM 


A 

A* 


2  3 

(1  +  .  2M  ) 
1.728M 


(This  relation  is  illustrated 
in  Figure  13) 


A  correction  is  then  applied  to  account  for  the  fact  that  the 
effective  throat  area  is  smaller  than  the  geometric  throat 
area.  This  correction  is  presented  in  the  form  shown  in 
Figure  14.  The  correction  is  based  on  correlations  of 
experimental  data  with  theory  used  as  a  guide  in  fairing 
the  curves  where  data  points  are  not  available.  If  the  M 
is  1.00,  then  the  assumed  recovery  is  reduced  by  .05  iEFF 
and  the  process  is  reiterated,  beginning  with  the  calculation 
of  A-/Ac  until  a  value  of  is  obtained  which  is  just  at 

or  below  M  =1.0.  iEFF 

EFF 


The  diffuser  loss  coefficient  is  then  obtained  from  input  data 
tables  as  a  function  of  M,r  .  An  example  of  this  input  data 

EFF 

is  presented  in  Figure  16. 

Using  the  diffuser  loss  coefficient  obtained  above,  the 
diffuser  total  pressure  loss  is  calculated  from: 


If  there  is  also  a  loss  in  total  pressure  from  froe-stream 
to  inlet  location,  the  two  losses  are  combined  to  yield 
the  final  total  pressure  recovery: 


If  the  calculated  recovery  differs  by  more  than  .001  from  the 
assumed  value,  then  the  process  is  reiterated  using  a  new  value 
of  recovery  which  is  the  average  of  the  calculated  and  assumed 
values . 


4.2.3  Supersonic 

The  major  causes  of  inlet  total  pressure  losses  in  the 
supersonic  speed  regime  (1.20  5  M  i  4.5)  are  shock  losses 
and  subsonic  diffuser  losses.  Shock-boundary  layer  interactions 
and  the  amount  and  distribution  of  boundary  layer  control 


43 


(1) 

Ref:  PG  117 

NA  68-655 

VoIX 

(2) 

Subton  ic  Entrance 

Profile 

(3) 

Diffuter  L/ht  r  9 

(4) 

Thin  Boundary  Layer 

Figure  16:  DIFFUSER  LOSS  COEFFICIENT 
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bleed  airflow  also  affect  both  shock  losses  and  subsonic 
diffuser  efficiency.  However,  it  is  difficult  to  separate 
out  the  contribution  of  boundary  layer  effects  on  either 
shock  losses  or  diffuser  efficiency  without  extensive  use 
of  experimental  data,  which  in  many  cases  does  not  exist. 

The  calculation  method  in  the  PITAP  procedure  uses  estimated 
inviscid  shock  losses  and  subsonic  diffuser  losses  based  as 
much  as  possible  on  the  results  of  experimental  work.  The 
effects  of  boundary  layer  on  shock  losses  are  accounted  for 
by  correcting  the  ideal  shock  losses  to  make  them  agree 
closer  with  experimental  results.  The  effect  of  boundary 
layer  on  subsonic  diffuser  losses  is  accounted  for  by  using 
realistic  diffuser  duct  loss  coefficients  obtained  from 
experimental  data. 

4. 2. 3.1  Two-Dimensional  Inlets 


Inviscid  total  pressure  recovery  of  an  inlet  supersonic 
diffuser  depends  on  the  number  and  strength  of  shocks  used 
to  decelerate  the  flow  from  supersonic  to  subsonic  Mach  numbers. 
Figure  17  shows  the  trends  in  total  pressure  recovery  as  a 
function  of  free-stream  Mach  number  and  the  number  of  shocks, 
including  all  oblique  shocks  and  the  normal  shock.  Also  shown 
in  Figure  17  is  a  dashed  line  representing  the  approximate 
maximum  obtainable  recovery  based  on  experimental  data.  When 
an  inlet  is  operating  at  a  free-stream  Mach  number  which  is 
at  or  below  its  design  Mach  number,  and  has  variable  ramps 
that  can  be  positioned  to  provide  nearly  optimum  ramp  geometry, 
the  total  pressure  recovery  can  be  obtained  from  the  curves 
of  Figure  17.  The  maximum  practical  recovery  (shown  as  a 
dashed  line  in  Figure  17)  is  used  as  an  upper  bound.  For 
the  case  where  an  inlet  is  operating  above  its  maximum  design 
Mach  number  (except  for  a  normal  shock  inlet,  N  =  1) ,  flow 
stability  problems  are  commonly  encountered  due  to  shock 
waves  entering  the  inlet.  Therefore,  this  condition  is  not 
of  much  practical  interest  for  most  preliminary  design  studies. 

The  ramp  geometry  required  to  provide  the  estimated  total 
pressure  recovery  can  be  obtained  from  Figure  18,  which  shows 
the  total  deflection  for  optimum  shock  systems  having  various 
numbers  of  oblique  shocks  as  a  function  of  free-stream  Mach 
number.  If  all  the  compression  were  accomplished  externally, 
the  summation  of  the  deflection  angles  would  be  equal  to  the 
last  ramp  angle.  If  only  a  portion  of  the  compression  is 
done  externally,  the  last  external  ramp  angle  is  correspondingly 
less.  In  this  way,  depending  on  the  specified  external/internal 
compression  split,  the  last  external  ramp  angle  can  be  calculated 
from  the  equation: 
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Figure  18:.  Shock  Flow  Deflections 
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where 

split. 


(l  6  ex  A  . 

r*» — 7  1 


is  the  specified  external/total  compression 


The  inlet  throat  size  is  calculated  from  the  following  equation: 


(VAc>  .  /a\  /a*t> 


*Ao^Ac  ^ 


(9) 


where 


=  throat  area,  in. 


=  inlet  capture  area,  in. 


Ao/Ac 


A/A* 

A$/A5 


inlet  mass  flow  ratio  passing  through 
the  throat 

(1  +  *iMT2)3 


T775TH 


t 


=  l 


*x  /pt 


A  flow  chart  showing  the  steps  involved  in  the  calculation 
procedure  is  presented  in  Figure  62. 

4. 2. 3. 2  Axisymmetric  Inlets 

The  total  pressure  recovery  of  an  axisymmetric  inlet  can  be 
obtained  from  the  plotted  data  of  Figures  19  and  20  in  a 
similar  manner  to  that  used  for  two-dimensional  inlets.  Figure 
19  presents  recovery  plots  for  external  compression  inlets 
and  Figure  20  presents  plots  for  inlets  with  external  and 
internal  compression. 

4.3  BOUNDARY  LAYER  BLEED  DRAG 

To  calculate  the  drag  of  boundary  layer  bleed  air,  it  is 
necessary  to  know  the  geometry  and  location  of  the  bleed 
exit  nozzle,  the  amount  of  bleed  airflow,  and  the  total 
pressure  recovery  of  the  bleed  air  at  the  nozzle  exit. 

Methods  for  estimating  the  last  two  parameters  and  calculating 
the  drag  are  described  in  this  section. 
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Figure  20:  MAXIMUM  INV/SCID  TOTAL  PRESSURE  RECOVERY  FOR  AXISYMME  TRIC 
INLETS  ( EXTERNAL  AND  INTERNAL  COMPRESSION ) 


4.3.1  Airflow 


The  amount  of  bleed  air  required  is  obtained  from  Figure  21. 
These  data  were  obtained  from  experimental  data  (References  2 
through  6) ,  and  represent  a  nearly  optimum  bleed  configuration. 

4.3.2  Recovery 

The  bleed  plenum  total  pressure  recovery  is  obtained  from 
Figure  22.  These  pressure  recoveries  represent  average 
values  obtained  from  experimental  data  (References  3  through 
5)  . 


A  flow  chart  showing  the  steps  in  the  procedure  required  to 
obtain  bleed  airflow  and  total  pressure  recovery  is  presented 
in  Figure  63. 


4.3.3  Drag 

The  output  of  the  boundary  layer  bleed  subroutine  is  used  as 
input  to  the  momentum  and  flap  drag  subroutines  which  are 
used  for  both  bleed  drag  and  bypass  drag  calculations.  The 
only  difference  is  that  the  bleed  air,  A  /A  ,  is  used  in 

BLC 

bleed  calculations  and  bypass  air  is  used  in  bypass  drag 
calculations. 

Bleed  (and  bypass)  drag  is  composed  of  two  parts:  momentum 
drag  and  flap  (or  exit)  drag.  Momentum  drag  depends  on  the 
exit  angle  of  the  discharged  air,  total  pressure  recovery  of 
the  discharged  air,  type  of  nozzle  through  which  the  air  is 
exited,  and  free-stream  Mach  number.  Flap  drag  depends  on 
door  angle,  door  aspect  ratio,  and  free-stream  Mach  number. 


4. 3. 3.1  Momentum  Drag 

The  procedure  for  calculating  momentum  drag  is  presented  in 
the  flow  chart  of  Figure  64.  The  input  quantities  required 
for  the  calculation  are: 


local  Mach  number,  dimensionless. 

exit  angle  of  air  relative  to  free-stream, 
degrees . 

total  pressure  recovery  of  bypassed  airflow, 
dimensionless . 


boundary  layer  bleed  mass  flow  ratio, 
dimensionless . 
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Figure  22  PL  ENUM  TO  TA L  PRESSURE  RECO  VER  Y 


bypass  mass  flow  ratio,  dimensionless. 


2 

A^  ~  inlet  capture  area,  in. 

Nozzle  type  must  be  specified  as  one  of  the  following  types: 

(a)  ~  Convergent  nozzle 

(b)  ~  Convergent-divergent  nozzle,  fully  expanded 

(c)  ~  Convergent-divergent  nozzle  with  a  specifier 

area  ratio,  A^/A^, ,  dimensionless 


A„  =  nozzle  exit  area,  in.2 

2 

At  =  nozzle  throat  area,  in. 

(a)  Convergent  Exit 

The  first  step  in  the  calculation  of  momentum  drag  for  the 
convergent  nozzle  is  a  test  to  determine  if  the  nozzle  is 
choked.  This  is  accomplished  by  calculating  M0  according 
to  the  following  equation:  min 

(10) 


Mq  represents  the  lowest  free-stream  Mach  number  at 
min 

which  the  available  pressure  ratio  will  be  sufficient  to 
choke  the  nozzle  throat.  If  MQ  £  MQ  ,  the  nozzle  throat  is 

assumed  choked,  Mg  =  1.0,  and  Pg/P0  =  1.0. 

Next,  the  nozzle  throat  area  is  calculated  by  the  equation: 
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Drag,  as  a  fraction  of  the  free-stream  momentum  of  the  bypassed 
airflow,  D  ,  is  calculated  by  the  following  equation: 


(WaV 


(12) 
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Jj 

Converting  the  above  drag  to  CD: 
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(13) 


If  exit  area,  A„,  and  P-,  /P,_  are  specified  instead  of  A„  /h„ 
E  TE  T0  °BLC  C 

and  P  /P  ,  (as  above)  the  following  calculation  procedure  is 
lE  0 

used  to  determine  bleed  drag: 

The  bleed  airflow,  A  /A  ,  is  first  computed  (if  it 

BLC  C 

is  not  specified)  by  the  following  equation: 
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BLC 
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p 
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After  An  /Ap  is  calculated,  the  drag  calculation  proceeds  in 
BLC 

the  same  way  as  before,  using  Equation  (12). 


For  an  unchoked  convergent  nozzle,  (M 


0  Mo 


) ,  the  exit 


nun 


static-to-total  pressure  ratio,  P../P  ,  is  calculated  ly  the 

E 

following  steps: 

P~  _ l _ 


~3T5 
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(1  +  .2  Mq  ) 
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Next,  the  exit  static-to- total  pressure  ratio  is  calculated: 


(15) 


Next,  M„  is  obtained  from  P_  using  Equation  (14). 
E  E 


A 

Then  (rrj-)  is  obtained  using  Mp  and  the  following  equation: 


(A-) 

'a*'. 


(1  +  .2  ) 

1.728  M~ 


The  throat  area  is  then  calculated  by: 


Drag,  as  a  fraction  of  free-stream  momentum,  is  calculated 
by: 


(17) 
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a  0 


=  1  -  cos  9 
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CD  is  then  calculated  as  before: 
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(b)  Convergent-Divergent  Exit 

A  test  is  first  performed  to  determine  whether  the  nozzle 
throat  is  choked.  MQ  is  calculated  and  compared  with  MQ 

min 

as  in  the  case  of  the  convergent  noz?le  calculation.  If  the 
calculation  shows  that  the  nozzle  throat  is  not  choked,  the 
C-D  nozzle  is  not  a  good  choice  and  the  configuration  should 
be  re-examined,  since  excessive  drag  would  result  from  the 
use  of  this  configuration.  If  this  happens,  the  computer 
program  terminates  the  case  and  goes  to  the  next  case. 

If  the  test  shows  that  the  throat  is  choked,  two  calculation 
procedures  can  be  followed. 


(c)  C-D  Exit  Fully  Expanded 


If  it  is  specified  in  the  input  ita  that  the  C-D  nozzle  is 
to  be  fully  expanded,  the  exit  static  pressure,  PE,  is  assumed 

to  be  equal  to  ambient  pressure,  PQ,  and  the  nozzle  exit  static- 

to-total  pressure  ratio  P  /P  ,  is  calculated  from  the  following 

LE 

equation: 


Using  the  above  ratio,  exit  Mach  number  is  obtained  from  the 
relation: 


1 

(1  +  .2Me2) 


3.5 


Next,  the  free-stream  static-to-total  pressure  ratio  is 
calculated: 


3.5 


(1  +.2  Mq  ) 
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Then  P„  is  calculated  from  the  relation: 


Nozzle  area  ratio  is  calculated  by: 

Ae  (1  +  .2  Me2)  3 

=  1.728  Me 

Nozzle  throat  area  size  is  calculated  by: 


Nozzle  exit  area  is  calculated  from: 


(18) 


(19) 


(20) 


Drag  is  then  computed  from  the  following  relation: 


(21) 


58 


Converting  to  CD: 


(d)  Convergent-Divergent  Exit  with  Specified  Area  Ratio 

If  the  throat  of  the  nozzle  is  choked,  as  will  have  been 
determined  from  the  test  that  was  made  using  MQ  >  MQ  ,  the 
calculation  proceeds  as  described  below.  min 

The  nozzle  area  A£/A*  is  set  equal  to  the  specified  nozzle  area 

ratio  AE/AT.  Then  the  .^ree-stream  (A/A*)Q  is  obtained  from 

Figure  13.  These  data  provide  the  necessary  items  to  calculate 
nozzle  throat  area  according  to  the  following  equation: 


Then  nozzle  exit  area  is  computed  from: 


AE 


\ 

The  nozzle  exit  Mach  number,  M£,  is  obtained  from  Figure  13 
as  a  function  of  the  nozzle  A£/A*  ratio.  Then  the  exit 
pressure  ratio,  P£/PT  >  is  calculated  from  the  equation: 
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Also,  the  free-stream  pressure  ratio  is  calculated  from  the 
same  equation. 


Next,  the  exit  pressure  ratio,  PE 


is  calculated  from  : 


Drag  is  then  computed  from: 


This  is  converted  to  drag  coefficient  by: 


All  the  equations  for  bypass  momentum  drag  and  flap  drag  are 
the  same  as  those  used  for  boundary  layer  control  bleed  airflow 
drag  calculations. 

4. 3. 3. 2  Flap  Drag 

Flap  drag  is  calculated  by  integrating  the  predicted  pressure 
coefficients  acting  on  the  exit  doors  in  the  flight  direction. 

To  accomplish  this,  tables  of  average  pressure  coefficients  for 
doors  of  various  aspect  ratios  (W/h)  from  .75  to  4.0  are  used  as 
input  data  to  obtain  the  pressure  coefficients. 


£0 


A  correction  is  also  applied  to  account  for  the  fact  that  the 
doors  are  usually  operating  in  a  region  of  turbulent  boundary 
layer.  The  detailed  steps  used  in  the  calculation  procedure 
are  shown  in  Figure  64  and  described  in  the  sections  which 
follow. 


w 

door  width,  inches 

h 

length  of  door,  inches 

n 

number  of  doors,  dimensionless 

A0  /A 

Bleed  mass  flow  ratio 

ac 

,  .  2 

capture  area  size,  m 

Mo 

free-stream  Mach  number,  dimensionless 

PT 

ae 

PT 

0 

total  pressure  recovery  of  the  airflow, 
dimensionless 

cp  vs  M0' 

&  0  eff 

tables  of  average  pressure  coefficient 

free-stream  Mach  number,  door  aspect  ratio, 
and  door  angle. 

Door  aspect  ratio  is  computed  from  the  input  door  dimensions: 

A<  -  W/h  (23) 

The  next  step  in  the  calculation  procedure  is  to  compute  the 
ratio  of  door  exit  area  to  capture  area,  A^/A^ .  Two  possibilities 

exist:  convergent  nozzle  and  convergent-divergent  nozzle. 

For  a  convergent  exit,  the  exit  area  is  calculated  by  the 
following  equation: 
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(where  A£/A*E  =  1.0  for  a  choked 
convergent  nozzle) 


(  °BLC\ 

Vac  1 


(.95) 

(0.95  is  the  assumed  flow 
coefficient) 


(24) 


This  equation  assumes  that  the  exit  consists  of  n  doors  of  equal 
size. 

The  door  flap  angle  is  calculated  from: 


(0FLAP) 


*  TAN 
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Where  aflap/ac 


Wh 


(26) 


Using  (8  PTap)  and  free-stream  Mach  number,  M  ,  the 

ACTUAL  0 

effective  door  angle  is  obtained  from  the  plotted  data  of 
Figure  23,  which  shows  the  effect  that  turbulent  boundary 
layer  has  on  reducing  flap  angle.  From  the  plotted  data, 
(0PTAp)  •  the  effective  door  flap  angle,  is  obtained  as  a 
EFF 

function  of  actual  door  angle  and  free-stream  Mach  number, 
assuming  average  turbulent  boundary  layer  conditions. 

The  average  pressure  coefficient  over  the  surface  of  the  door 
is  obtained  from  the  data  of  Figure  24  as  a  function  of 

(0FLAP)Epp'  M0'  and 

Using  the  pressure  coefficients,  the  flap  drag  coefficient  is 
calculated  from  the  following  equation: 
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(27) 


62 


piap)  *'»  Effective  Flap  Angle 
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v  24.  BYPASS  DOOR  PRESSURE  COEFFICIENTS  (Cont) 


Figure  24  B  YPASS  DOOR  PRESSURE  COEFFICIENTS  (Corn) 
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If  a  convergent-divergent  exit  nozzle  is  used,  the  nozzle 
throat  area  is  first  calculated  by: 


(Assumes  (A/A*),f  =  1.0 
for  a  choked  throat) 


77§5) 


(28) 


The  nozzle  exit  area  is  then  obtained  from  input  curves  of 

A 

E  vs  A_  which  represent  the  geometric  characteristics  of 

A  1 

aT 

the  doors.  Exit  area  is  then  calculated  from  the  relation: 


The  calculation  then  proceeds  the  same  way  as  for  a  convergent 
nozzle. 

4 . 4  BYPASS  DRAG 


The  calculation  procedure  for  computing  bypass  drag  is  the 
same  as  that  for  boundary  layer  bleed  drag  except  that  bypass 
total  pressure  recovery  is  different  (usually  higher)  than 
that  of  the  boundary  layer  bleed  air.  The  recovery  of  the 
bypass  air  may  be  assumed  to  be  .85  x  (P  /P  )  for  preliminary 

2  *0 

studies.  If  a  more  refined  calculation  is  required,  the  bypass 
recovery  may  be  varied  linearly  from  a  value  of  1.00  x  ( P  /P  ) 

with  the  bypass  doors  closed  to  a  value  of  .70  x  (P  /P  )  with 

l2  lO 

the  doors  full  open  (which  usually  bypasses  around  fifty  per  cent 
of  capturea  inlet  airflow) .  The  recommended  bypass  recovery 
factors  are  illustrated  in  Figure  25. 

The  method  described  in  previous  sections  for  calculating  exit 
drag  for  bypass  and  boundary  layer  bleed  air  has  covered  only 
the. specific  case  where  flap-type  doors  are  used  for  exits. 

The  same  methods  can  also  be  used  for  other  types  of  exits 
by  use  of  the  engineering  assumptions  shown  in  Table  III. 


68 


Figure  25:  BYPASS  AIRFLOW  TOTAL  PRESSURE  RECOVERY 


Table  III.  ENGINEERING  ASSUMPTIONS  FOR  BOUND  A  R  Y 
LA  YER  BLEED  AND  BYPASS  EXITS 


Type  of  Exit 


Assumption 


Use  Method 
Described 
in  Volume  I 


Slot-Type  Exits 


Assume  Exit 
Drag  is  Same 
as  for  Flap 

8t  aAir  “  aFlap 
Because  External 
Pressure  Field 
Will  be  Created 
Just  as  if  a  Dooi 
Were  Present 


4.5  SPILLAGE  DRAG 


Spillage  drag  is  calculated  from  the  following  equation: 


'D 


=  K 


SPILL 


ADD 


ADD 


It  is  defined  as  the  incremental  change  in  total  airplane 
drag  due  to  spilling  excess  air  ahead  of  the  inlet.  However, 
it  is  usually  measured  as  an  increment  (AC  )  from  the 

dPILL 

drag  level  obtained  for  the  full  airplane  configuration  with 
the  inlet  operating  at  a  baseline  mass  flow  ratio.  When  the 
inlet  is  operating  at  this  baseline  mass  flow  ratio,  the 
ACn  is  zero.  The  spillage  drag  at  the  baseline  mass 
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flow  ratio  is  included  in  the  airplane  aerodynamic  drag.  The 
baseline  mass  flow  ratio  is  specified  for  each  inlet  configuration 
as  a  function  of  f ree-stream  Mach  number.  It  is  normally 
selected  at  mass  flow  ratios  where  spillage  effects  will  be 
a  minimum,  for  realistic  operating  airflow  conditions.  This 
baseline  mass  flow  is  normally  taken  to  be  the  critical  mass 
flow  ratio.  At  supersonic  speeds,  with  the  shock  wave  attached 
to  the  first  ramp,  critical  inlet  operation  means  that  the 
normal  shock  is  at  or  inside  the  cowl  lip.  lor  subsonic  and 
detached  shock  wave  conditions,  it  means  that  inlet  captured 
free-stream  tube  area,  AQ,  is  equal  to  the  physical  flow  area 

at  the  cowl.  The  establishment  of  the  baseline  mass  flow 
ratio  as  described  provides  a  basis  for  accounting  of  aero  and 
propulsion  forces.  The  throttle-dependent  drag  is  thus  included 
in  the  spillage  drag  (which  is  accounted  for  in  installed  net 
thrust)  and  the  drags  that  are  independent  of  throttle  setting 
are  included  in  the  aerodynamic  drag  of  the  airplane. 

4.5.1  Theoretical  Additive  Drag 

CQ  is  the  theoretical  additive  drag  of  the  inlet.  It  is 
ADD  computed  by  several  different  methods,  depending  on 
the  configuration  of  the  inlet,  free-stream  Mach  number,  and 
shock  geometry.  This  section  describes  the  methods  which  are 
used  for  each  particular  situation. 


The  additive  drag  of  an  open-nose  inlet  is  calculated  as 
shown  in  Figure  26.  This  equation  can  be  used  for  both 
subsonic  Mach  numbers  and  supersonic  Mach  numbers  where  a 
normal  shock  is  standing  ahead  of  the  inlet. 

For  two-dimensional  inlets  at  subsonic  Mach  numbers  and  for 
supersonic  Mach  numbers  with  the  terminal  shock  at  or  inside 
the  inlet  lip,  the  method  shown  in  Figure  27,  is  used. 


Figure  26.  ADDITIVE  DRAG  OF  A  PITOT  INLET 


Figure  2 . ...  ADDITIVE  DRAG  OF  A  TWO/DIMENSIONAL  INLEi 


The  last  term  on  the  right  hand  side  of  the  equation  (ramp 
drag)  represents  the  incremental  change  in  ramp  drag  due  to 
reducing  the  mass  flow  ratio  below  the  reference  value.  This 
incremental  drag  is  obtained  by  calculating  the  pressure  drag 
on  the  ramp,  using  the  average  pressure  on  the  ramp,  (p0+P1)/2 

and  the  projected  ramp  area,  A^.  The  reference  ramp  drag  is 

obtained  from  the  plotted  data  of  Figure  28,  as  a  function  of 
ramp  thickness  ratio,  t/c,  and  free-stream  Mach  number. 

For  two-dimensional  inlets  operating  at  supersonic  Mach  numbers 
with  the  normal  shock  outside  the  inlet  lip  (subcritical) , 
the  additive  drag  calculation  method  of  Reference  7  is  used. 

This  method  has  been  programmed  for  the  CDC  6600  computer  as 
a  separate  program  and  is  not  included  in  this  report. 

Theoretical  additive  drag  for  an  axisymmetric  inlet  at  subsonic 
Mach  numbers  and  for  supersonic  Mach  number  conditions  where 
the  normal  shock  is  at  or  inside  the  inlet  lip,  the  method 
shown  in  Figure  29  is  used. 

Incremental  spike  drag  is  calculated  in  much  the  same  way  as 
ramp  drag  for  the  two-dimensional  inlet.  The  average  spike 
pressure  and  projected  spike  area  are  used  to  calculate  the 
incremental  drag.  The  spike  reference  drag  is  obtained  from 
pressure  drag  data  on  cones  (Figure  30) . 

For  axisymmetric  inlets  operating  with  the  terminal  shock 
ahead  of  the  inlet  lip,  the  method  documented  in  Reference  8 
is  used.  This  method  has  been  programmed  for  the  CDC.  6600 
computer  as  a  separate  program  and  is  not  documented  in  this 
report. 

4.5.2  K.„  Factors 

ADD 


The  theoretical  additive  drag  is  multiplied  by  a  correction 
factor,  Kadd,  based  on  experimental  data,  to  account  for  the 

configuration  effects.  These  effects  include  the  effect  of 
cowl  lip  shape,  bluntness,  and  sideplate  cutback.  The  K  ^ 
factors  are  obtained  by  comparing  experimentally 
measured  drag  variation  as  a  function  of  inlet  mass  flow  ratio 
with  the  theoretically  calculated  drag  variation.  A  catalog 
of  factors  to  be  used  is  presented  in  Volume  IV. 

4 . 6  ANGLE-OF-ATTACK  EFFECTS 

Correction  procedures  are  discussed  here  for  the  three  most 
frequently  encountered  types  of  inlet  installations:  free- 
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from  Ref.  7 


1(7  +  1)  t/c) 


Figure  29.  ADD/TI VE  DRAG  OF  AN  AX/SYMMETR/C  INLET 
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Figure  30:  PRESSURE  DRAG  ON  CONES  ,7V  TRANSONIC  FLOW 
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stream,  underwing  and  fuselage  mounted. 

For  free-stream  installations,  each  angle  of  attack  position 
is  calculated  as  a  separate  data  point  without  applying  a 
correction  for  the  presence  of  the  wing  or  the  forebody. 

For  all  types  of  installations,  it  is  assumed  that  the  basic 
inlet  is  installed  in  a  location  where  it  is  free  of  adverse 
boundary  layer  effects . 

For  the  underwing  inlet  location,  it  is  assumed  that  the 
inlet  is  located  within  the  flowfield  of  a  two-dimensional 
wing.  As  angle-of-attack  is  increased,  the  wing  compression 
field  provides  a  local  flowfield  at  the  inlet  that  is 
different  from  free-stream  according  to  the  relations  plotted 
in  Figures  31  through  34. 

Figure  31  presents  a  correction  factor,  K,  used  to  compute 
the  local  Mach  number  behind  a  wing  shock.  It  is  derived 
from  the  linear  portions  of  the  vs  a  data  shown  in  Figure 
32.  u 

Figure  33  shows  the  effect  of  angle-of-attack  on  local  stream 
tube  area  for  supersonic  Mach  numbers. 

For  preliminary  design  purposes,  the  performance  calculations 
for  an  inlet  in  the  wing  flow  field  are  approximately  equal 
to  those  for  an  inlet  in  free-stream  except  that  the  effect 
of  wing  compression  on  local  inlet  flowfield  properties  must 
be  accounted  for  before  the  inlet  performance  calculations 
are  made.  Figures  31-34  are  used  for  this  purpose. 

For  side-mounted  inlet  installations,  the  available  experi¬ 
mental  data  (Volume  IV)  indicates  there  is  relatively  little 
effect  of  angle-of-attack  on  local  inlet  Mach  number  and 
total  pressure  recovery  up  to  20  degrees  angle-of-attack. 
Therefore,  for  preliminary  studies,  it  is  assumed  that  the 
inlet  local  flowfield  properties  are  the  same  as  free-stream. 

For  under-fuselage  mounted  inlets,  the  fuselage  provides 
increasing  compression  as  angle-of-attack  increases.  This 
compression  field  reduces  the  local  Mach  number  ahead  of  the 
inlet  cind  can  also  reduce  the  total  pressure  recovery  ahead 
of  the  inlet. 

Although  this  compression  field  is  dependent  on  the  shape  of 
the  forebody  ahead  of  the  inlet  and  the  location  of  the  inlet, 
aft  of  the  forebody  nose,  the  experimental  data  have  prevailing 
trends  that  indicate  the  flowfield  effects  are  generally 
somewhere  between  those  produced  by  a  conical  forobody  and  a 
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Fig ,n  33:  EFFECT  OF  ANGL E  OF  A  TTACK  ON  STREAM  TUBE 
6 REA  FORA  TWO-DIMENSIONAL  WING 
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two-dimensional  surface.  The  correction  curves  to  be  used 
for  the  under-fuselage  inlets  are  presented  in  Figures  35 
and  36.  They  are  based  on  experimental  data  from  Reference  9. 


Figure  36:  EFFECT  OF  ANGLE  OF  ATTACK  ON  LOCAL 
INLET  TOTAL  PRESSURE  RECOVERY  FOR 
AN  UNDER  FUSELAGE  MOUNTED  INLET 
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SECTION  V 


NOZZLE /AFTERBODY  PERFORMANCE 


5 .1  NOZZLE  INTERNAL  PERFORMANCE 

Figure  37  presents  a  diagram  showing  how  the  nozzle/ 
afterbody  calculation  procedures  relate  to  the  installed 
propulsion  system  performance  procedures.  Each  of  the 
nozzle/afterbody  procedures  is  discussed  in  detail 
in  the  following  sections. 

5.1.1  Convergent 

Test  data  are  available  in  Volume  IV,  Section  III  that 
show  the  measured  velocity  coefficients  (Cy)  for  a  series 
of  convergent  nozzles  as  a  function  of  entrance  to  exit 
diameter  ratio,  pressure  ratio,  and  nozzle  convergence 
angle.  The  velocity  coefficients  are  defined  as  follows: 

Cv  »  (FG) MEASURED _ 

^ ACTUAL  (V) ISENTROPIC,  FULLY  EXPANDED 

These  measured  C„s  are  used  with  the  ideal,  fully 
expanded  gross  truust  calculated  by  methods  described 
in  Section  III  to  obtain  the  predicted  nozzle  gross 
thrust  for  a  given  mass  flow.  Convergent  nozzle  dis¬ 
charge  coefficients  were  also  measured  during  the  tests. 
These  coefficients  presented  in  Reference  29, are  used  to 
relate  the  effective  to  actual  throat  areas  for  purposes 
of  sizing  the  nozzle  throat  to  pass  a  given  mass  flow  or, 
conversely,  for  determining  the  effective  throat  area  for 
a  given  throat  geometric  area. 

The  experimental  C„s  are  accurate  to  within  +  .25%,  and 
should  be  adequateror  most  preliminary  design  purposes. 

If  it  is  desired  to  calculate  the  internal  performance  of 
a  convergent  nozzle  by  theoretical  buildup,  it  is 
recommended  that  Cy  be  calculated  as  follows: 
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INLET  PERFORMANCE  INSTALLED  PROPULSION  |  NOZZLE/AFTERBODY 

CALCULATION  SYSTEM  PERFORMANCE  i  CALCULATION 


NOZZL E/A F TERBOD Y  CAL CULA  TiON  PROCEDURE 


s 


Cy  -  I* 


AC, 


V, 


1.  - 


W8  +A8  «VW 


where  Ag  = 


= 


actual  throat  area  =  Ai/Cn  ;  Ai  =  thermodynamic 

4  throat  area, 

in  .  , 

experimental  throat  discharge  coefficient, 
dimensionless  (Ref.  43) 

2 


Pg  =  static  pressure  at  the  throat,  lb/in. 


Vg  =  jet  velocity  at  the  throat,  ft/sec 


ideal  gross  thrust  for  fully  expanded  flow 

=  internal  loss  due  to  skin  friction  (cal¬ 
culated  by  methods  of  Ref.  44) 


Theoretical  methods  are  also  contained  in  Reference  44  for 
calculating  the  effect  of  boundary  layer  growth  on  nozzle 
discharge  coefficient. 

5.1.2  Convergent-Divergent 

The  performance  of  a  convergent-divergent  nozzle  is  determined 
from  the  stream  thrust  coefficient  discussed  in  Volume  IV, 
Section  III. 


The  procedures  for  calculating  the  thrust  of  a  convergent- 
divergent  nozzle  are  given  in  Section  III.  By  these  procedures, 
the  gross  thrust  is: 


F 

-  cs  (j£ 

S  g 

v9  +  p9a9) 

-  p  ,  a„ 

g 

amb  9 

where  Cg 

=  stream 

thrust  coefficient  (l-ACc) 

=  f  (Ag,R9,9)! 

from  Figure 
38.  Dimen¬ 

A9 

=  nozzle 

exit  area, 

in2 

sionless 

e 

=  nozzle 

divergence 

angle,  degrees 

5.1.3  Ejector 

There  is  not  sufficient  test  performance  data  in  Volume  IV, 
Section  5  to  calculate  ejector  nozzle  performance  fcr  all 
operating  conditions.  The  discussion  in  Reference  10  indicates 
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Figure  38:  CONICAL  CONVERGENT/DIVERGENT  NOZZLE  NON-ISENTROPIC  EXPANSION  LOSS 


that  for  certain  conditions  ejector  nozzle  performance 
calculated  by  compound  flow  theory  and  experimental  results 
show  excellent  agreement.  Therefore,  this  method  was  selected 
for  use  in  calculating  ejector  performance  where  no  test 
data  are  available.  Using  this  method,  the  nozzle  thrust  for 
an  ejector  nozzle  with  dual  flow  can  be  calculated  as  follows: 


W  V 
S  9, 


+  P9A*. 


IT  +  P9A9, 


P 


ambA9 


where  CA  =  divergence  loss  coefficient  =1/2  (l+Cos  0) , 
w  dimensionless 

Wg  =  secondary  flow  rate,  3b/sec 
Wp  =  primary  flow  rate,  lb/ sec 

2 

Aq  a  nozzle  exit  area  for  the  primary  stream,  in. 

Aq  =  nozzle  exit  area  for  secondary  stream,  in.^ 

S 

.  2 

A(j  =  nozzle  exit  area,  m. 


V.  =  nozzle  exit 


velocity  of  secondary  airflow,  ft/sec 


Vq  =  nozzle  exit 
P 


velocity  of  primary  airflow, 


f t/sec 


The  nozzle  geometry  and  secondary  inlet  flow  conditions  must 
be  specified  to  calculate  the  secondary  flow  (w^) .  If  compound 
choking  does  not  occur  (see  Reference  10)  tnen  ‘  the  secondary 
flow  is  calculated  at  the  exit  plane  by  iterating  until  the 
static  pressures  of  beth  streams  are  matched  (See  Section  3.2. J. 3). 


P 


W 


s 


55  53.34 2  ' Tq 

y6 


The  sane  procedure  is  repeated  for  the  nozzle  throat  and  if 
the  calculated  secondary  flow  is  less  than  the  exit  flow,  it 
is  assumed  that  this  exit  area  can  be  reduced  to  mutch  the.  flow. 
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An  iteration  on  A^  to  obtain  maximum  secondary  is  required 
for  the  compound- choked  flow  condition. 

5.1.4  Plug 

There  is  not  sufficient  test  performance  data  in  Volume  IV, 

Section  III  to  calculate  plug  nozzle  performance  for  all 
operating  conditions.  A  good  estimate  of  plug  nozzle 
performance  can  be  made  by  subtracting  toe  added  skin 
friction  drag  due  to  the  external  plug  surface  from  the 
calculated  gross  thrust  of  a  concentric  nozzle: 


where 


(0.014) 


It  is  assumed  the 
to  fully  expanded 
pressure,  q,  over 
drag. 


dskin 


velocity  at  the  end  of  the  plug  is  equal 
isentropic  flow.  The  average  dynamic 
the  plug  is  used  to  calculate  skin  friction 


where  Cp  =  average  skin  friction  coefficient,  dimensionless. 

DSKIN  =  s^;'-n  friction  drag,  lb 
(4tt  Ag)1/2 

Dh  =  Perimeter  '  dimensionless 
e 

Cy  =  maximum  thrust  coefficient  for  a  concentric 
MAX  nozzle  (Reference  28)  dimensionless 

5.2  NQZZLE/AFTERBODY  DRAG 

Nozzle/afterbody  drag  includes  all  drag  elements  associated 
with  the  exhaust  system  installation  that  are  affected  by 
the  engine  power  setting.  Crag  associated  with  skin  friction 
is  included  in  the  airplane  drag  polars.  Drag  elements 
identified  as  contributors  to  the  nozzle/afterbody  drag 
include: 
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1.  Nozzle/afterbody  boattail  drag 

2.  Interference  effects  for  multi  engine  installations 

3.  Afterbody  base  drag 

4.  Scrubbing  drag 

Nozzle/afterbody  drag  is  computed  using  maps  which  represent 
the  afterbody  drag  characteristics  (Figure  39) ,  external  input 
geometric  data,  and  fundamental  engine  data  obtained  internally 
from  the  engine  subprogram.  The  external  inputs  are  required 
constants  that  describe  the  nozzle  and  afterbody  geometry: 
Maximum  nacelle  diameter,  boattail  length,  lateral  nozzle 
spacing,  base  area,  etc.  Fundamental  engine  data  obtained 
internally  from  the  engine  subprogram  include  nozzle  thrust 
area,  nozzle  pressure  ratio,  free-stream  conditions  and  ideal 
gross  thrust.  An  essential  geometry  input  is  the  nozzle 
exit  area,  A_,  which  is  required  for  boattail  drag  computation. 
This  parameter  is  obtained  xn  either  of  two  ways: 

1.  From  the  engine  subprogram  when  the  existing  nozzle 
data  are  used. 

2.  From  a  table  of  Ag/Ag  that  is  developed  when  new 
nozzle  data  are  used? 

The  logic  used  to  connect  the  input  drag  data  and  the  engine 
subprogram  is  illustrated  in  Figure  40. 

The  nozzle/afterbody  subprogram  requires  much  less  external 
input  data  than  the  inlet  subprogram.  This  is  due  to  the 
fact  that  the  maps  used  to  obtain  nozzle/afterbody  drags  have 
been  generalized  to  cover  a  wide  range  of  possible  configur¬ 
ations,  and  are  built  into  the  program.  Thus,  it  is  necessary 
to  supply  as  external  in. ut  data  only  the  geometric  parameters 
used  to  specify  the  afterbody  geometry. 

5.2.1  Boattail  Drag 

The  generalized  maps  which  are  now  in  the  program  are  shown 
in  Figures  41,  42,  43  and  44.  Figure  41  shows  the  variation 
of  boattail  drag  as  a  function  of  boattail  angle  and  Mach 
number  for  free-stream  Mach  numbers  up  to  .95.  From  Mach 
1.0  to  3.0,  the  equation  shown  in  Figure  41  is  used.  These 
drag  data  are  based  on  islatcd  wind  tunnel  data  for  circular 
arc  boattails  at  nozzle  pressure  ratios,  PT/P  of  2.5.  For 
nozzle  pressure  ratios  greater  than  2.5,  an  additional 
correction  is  applied  to  the  basic  boattail  drag  coefficient 
obtained  from  Figure  42.  This  added  correction  shown  in 
Figure  42  is  a  function  of  nozzle  pressure  ratio  and  boattail 
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Boattail 

Drag 


Interference 

Effect 


Derived  from  Isolated  Nozzle 
Wind  Tunnel  Test  Data 


Empirical  Correlation  of 
Lateral  Spacing  Effects  on 
Nozzle/Afterbody  Drag  Using 
Wind  Tunnel  Data 


Nozzle  Spacing  S/Dg 


Derived  from  Test  Data, 
Included  When  Base  Areas 
are  Present  in  the  Afterbody 
that  May  Be  Affected  by 
Engine  Power  Setting 


Scrubbing 

Drag 


Derived  from  Theory  and 
Empirical  Correlations 
Included  When  Nozzle 
Exhaust  Flow  Impinges  on 
Aircraft  or  Nacelle  Surfaces 


Figure  39:  NOZZLE/AFTERBODY  DRAG  REPRESENTATION 
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Internal  Inputs 


External  Inputs 


Figure  40:  NOZZLE/AFTERBODY  DRAG  SUBPROGRAM 
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Figure  4/ NOZZLE  BOA  TTAIL  PRESSURE  DRAG  COEFFICIENTS  AS 


Data  Sources:  1.  NASATMX-1960 

2.  Boeing  Test  Data 

3.  Unpublished  Boeing  SST  Data 
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Mach  Number 

Figure 43:  SUPERSONIC BOATTA/L  DRAG  CORRELATION 
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angle.  This  added  correction  procedure  is  included  in  existing 
nozzle/afterbody  subprogram. 

For  supersonic  Mach  numbers  greater  than  Mach  1.0,  PITAP  uses 
the  expression  given  in  Figure  41.  This  represents  an 
empirical  correlation  of  data  from  several  sources,  as 
shown  on  Figure  43.  The  correlation  was  biased  at  the  lower 
Mach  numbers  toward  data  obtained  with  blown  jet  simulation. 

The  form  of  the  correlation  was  suggested  by  supersonic  wave 
drag  considerations  (Reference  11)  and,  as  shown  in  the 
Figure  43,  it  does  a  reasonably  good  job  of  collapsing  the 
data. 

The  option  is  available  in  the  computer  program  to  calculate 
complete  afterbody  drag  as  a  function  of  the  ratio  of  nozzle 
exit  area  to  maximum  fuselage  cross-sectional  area,  A„/A.q. 

The  drag  maps  built  into  the  program  for  this  calculation 
are  shown  in  Figure  44. 

The  drag  of  a  plug  nozzle  subjected  to  subsonic  external 
flow  may  be  estimated  by  considering  the  cowl  and  the  external 
boundary  of  the  jet  to  define  an  equivalent  boattail  as  shown 
in  Figure  45.  The  jet  diameter  at  the  end  of  the  plug  should 
be  computed  assuming  full  expansion  to  ambient  static  pressure. 
The  drag  of  this  equivalent  boattail  may  then  be  estimated 
from  the  data  of  Figure  41.  Reference  30  contains  experimental 
data  comparing  measured  plug  nozzle  drag  with  C-D  nozzle  drags. 
This  data  substantiates  that  the  equivalent  nozzle  concept  is 
suitably  accurate  for  preliminary  design  analyses. 

The  drag  of  a  plug  nozzle  at  supersonic  speeds  and  high 
nozzle  pressure  ratios  is  simply  the  drag  of  the  solid 
boattail  created  by  the  nozzle  shroud.  This  drag  can  be 
estimated  from  the  data  of  Figure  41. 

The  installed  thrust-minus-drag  of  a  plug  nozzle  may  be 
estimated  by  summing  all  the  internal  and  external  losses 
for  the  nozzle  and  subtracting  these  losses  from  one.  Thus, 


F  -  D 
F. 

l 


D_ 

F  . 


=  1 


EACV 


'D 


^MAX  q«  *MAX 


IDEAL 


99 


Plug  Nozzle 


I 


Equivalent  Convergent  Nozzle 


Figure  45:  DEFINITION  OF  AN  EQUIVALENT  NOZZLE 


The  losses  in  internal  thrust  which  are  accounted  for  in  the 
velocity  coefficient  term,  IA C  ,  include  skin  friction,  non- 
ideal  expansion,  and  convergence  losses,  c  accounts 

for  boattail  drag.  Uamax 

5.2.2  Interference  Drag 

PITAP  is  programmed  to  accept  an  input  drag  increment  table 
giving  "interference"  drag.  Potentially  this  slot  could  be 
used  to  account  for  any  discrepancies  between  the  calculated 
isolated  nacelle  drag  and  the  actual  drag  of  the  nacelle 
installed  in  an  airplane  configuration. 

At  the  present  time  this  table  is  being  used  solely  to  account 
for  the  twinning  penalties  associated  with  putting  two  nacelles 
side  by  side  with  a  clean,  sharp-edged  interfairing  in  between. 
The  independent  variables  are  Mach  number  and  spacing  ratio, 
i.e.,  center  line  spacing  over  nozzle  exit  diameter.  (The 
present  version  of  PITAP,  however,  can  handle  any  number  of 
side-by-side  nozzles  in  a  row,  adding  n-1  twinning  penalties 
for  n  nozzles . ) 

The  interference  drag  coefficients  used  in  the  table  (Figure 
46)  at  spacing  ratios  greater  than  those  tor  which  C  is  a 
maximum  are  based  on  data  from  References  12  and  13,  I  which 
are  in  generally  good  agreement  with  recent  data  reported  in 
Reference  14.  At  lower  spacing  ratios,  the  curves  correspond 
to  completely  separated  flow  on  the  inside  half  of  the  nozzles 
and  on  the  entire  interfairing  and  base  drag  has  been  computed 
for  this  region. 

Figure  46  has  been  shown  to  work  quite  well  when  applied  to 
plug  nozzles  as  well  as  convergent  or  convergent-divergent 
nozzles.  In  this  case  the  equivalent  nozzle  concept  is 
first  used  to  estimate  isolated  nozzle  drag,  and  then 
Figure  46  is  used,  with  spacing  ratio  defined  to  be  center- 
line  spacing  over  fully  expanded  exit  diameters.  Figure  47 
(from  Reference  18)  shows  a  comparison  of  twin-nacelle 
thrust-minus-drag  predicted  performance  and  measured 
performance  for  both  plug  and  convergent-divergent  nozzles. 
Although  the  agreement  is  better  for  the  convergent-divergent 
nozzles,  in  both  cases  the  prediction  is  within  about  one  per 
cent  of  gross  thrust. 

For  nozzle  interfairings  that  differ  from  clean  sharp-edged 
fairings,  a  considerable  amount  of  data  (References  12,  14, 

15,  16  and  17)  is  available,  but  reliable  correlations  of  the 
interference  effects  as  such  do  not  exist. 
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Figure  46:  NOZZLE  INTERFERENCE  DRAG  COEFFICIENT 


Drag)/<deal 


In  principle  the  interference  table  now  in  PITAP  could  be 
replaced  by  an  additional  table  that  would  account  for  the 
interference  between  a  strut-mounted  engine  and  a  wing,  or 
perhaps  an  aft  fuselage.  A  generalized  table,  however, 
applicable  to  a  wide  range  of  configurations  would  be  extremely 
difficult  to  produce. 

In  a  recent  program  at  NASA  Lewis  (References  19,  20  and  21), 
interference  effects  between  an  underwing  engine  nacelle 
and  an  F106B  wing  were  evaluated  both  in  a  wind  tunnel  test 
and  in  actual  flight.  A  ,similar  test  with  a  simulated  wing 
is  also  reported  in  Reference  22.  These  reports  give  pressure- 
integrated  boattail  drags,  as  well  as  some  wing  pressure  data. 
Integrated  forces  on  the  wing  must  also  be  computed  before 
total  installation  loss  effects  and  interference  effects  can 
be  extracted  from  the  data.  These  tests  also  have  limited 
applicability  to  the  present  present  program  because  the 
emphasis  was  on  the  interference  effects  for  a  limited  range 
of  nozzle  geometric  variations  (e.g.,  only  a  15  degree  boattail 
angle  was  used)  rather  than  the  installation  losses  for  a 
range  of  simulated  power  settings. 

The  pressure  ratio  correction  used  in  PITAP  is  shown  in 
Figure  42.  The  effect  of  increasing  jet  pressure  ratio  is 
two-fold.  Increasing  the  size  of  the  exhaust  plume  tends 
to  force  higher  pressures  on  the  nozzle  surfaces,  decreasing 
drag.  Jet  entrainment  effects  tend  to  pump  down  the  pressure 
on  nozzle  surfaces,  particularly  where  separated  regions  exist, 
increasing  drag.  Over  a  moderate  range  of  pressure  ratios 
these  effects  tend  to  cancel  out.  Reference  23  provides  an 
example  where  these  two  opposing  effects  were  isolated.  Between 
jet-off  conditions  and  a  particular  jet  pressure  ratio  no  net 
effect  was  observed.  Using  a  solid  plume  simulator,  than 
eliminating  entrainment  effects,  it  was  determined  :  *  a  the 
jet  plume  produced  a  favorable  increase  in  isolated  o  ■  it  tail 
drag  coefficient  of  0.02,  while  entrainment  had  uncounted 
for  -0.02. 

The  relative  importance  of  entrainment  and  plume  effects  is 
a  strong  function  of  how  much  of  the  flow  is  separated  or 
in  a  base  region,  and  thus  sensitive  to  entrainment  effects. 

The  Boeing  data  of  Figure  48  clearly  illustrates  this 
phenomenon.  The  upper  three  bands  of  data  represent  twin 
nozzle  configurations  for  which  oil  flow  photographs  revealed 
significant  separated  regions.  Between  a  pressure  ratio  of 
two  (barely  more  than  ram  pressure  ratio)  and  four  there  is 
little  apparent  effect  of  pressure  ratio  on  drag.  The  dashed 
line  is  a  single  nozzle  reference  configuration  which  had  no 
separated  flow.  For  this  nozzle,  plume  effects  cause  a 
continuous  reduction  in  drag  with  increasing  pressure  ratio. 
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Similarly,  in  Figure  48,  the  highest  angle  nozzle  does  not 
benefit  as  much  from  high  pressure  ratio  as  the  intermediate 
angle  nozzles.  (The  low-angled  nozzles  are  almost  zero  drag 
at  the  outset.) 

Considerations  such  as  these  suggest  that  the  pressure  ratio 
effects  might  correlate  well  with  the  initial  drag  level  at 
a  reference  pressure  ratio.  If  the  drag  coefficient  were 
normalized  by  the  projected  nozzle/afterbody  external  surface 
area,  it  represents  a  measure  of  how  much  separation  existed 
at  the  reference  nozzle  pressure  ratio. 

5.2.3  Base  Drag 

The  PITAP  method  is  currently  set  up  to  receive  an  input 
table  of  base  pressure  coefficient  versus  Mach  number  for  a 
specified  afterbody  geometry.  It  can  use  the  base  pressure 
coefficient  from  the  input  table  and  a  specified  base  area 
to  compute  base  drag  as  a  function  only  of  free-stream  Mach 
number.  However,  to  be  throttle-sensitive,  the  base  area 
must  be  located  where  it  is  also  affected  by  the  propulsive 
jet  effects  which  vary  with  nozzle  pressure  ratio.  In  addition, 
the  base  pressure  is  known  to  be  sensitive  to  approaching 
boundary  layer  conditions.  Therefore,  methods  to  account 
for  these  effects  should  be  incorporated  into  the  calculation 
procedure  to  obtain  an  improved  base  pressure  prediction. 

Subso.nically,  the  effects  of  the  propulsive  jet  may  be 
accounted  for  by  using  correlations  obtained  by  McDonald 
and  Hughes  (Reference  24). 

The  McDonald  and  Hughes  correlation  gives  the  increment  in 
base  pressure  coefficient  between  a  base  with  no  jet  effects 
and  one  with  a  jet  of  a  given  nozzle  pressure  ratio  and  a 
given  temperature.  The  correlation  is  presented  in  Reference 
24  for  a  free-stream  Mach  number  of  0.9  and  jet  total  pressure 
ratios  of  2.0  and  3.0.  The  correlation  is  seen  to  work  for 
approach  boattail  angles  from  8  to  24  degrees,  and  specifically 
treats  the  case  of  an  annular  base  around  the  jet.  The  key 
correlating  parameter  is  the  jet  diameter  squared  divided  by 
the  product  of  the  base  diameter  and  the  maximum  diameter. 

Figure  49  presents  this  correlation  for  a  jet  pressure  ratio 
of  2.0.  The  success  of  this  parameter  in  correlating  the 
data  at  Mach  0.9  suggests  that  it  be  tried  at  other  subsonic 
Mach  numbers.  The  method  of  McDonald  and  Hughes  for  Mach  0.90 
is  now  programmed  into  PITAP. 

For  supersonic  Mach  numbers,  from  1.0  to  2.2,  an  excellent 
base  pressure  correlation  has  been  reported  in  Reference  25. 
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Tj  =  Propulsive  Jet  Total  Pressure 
Poo  =  Freest  ream  Static  Pressure 


Figure  49:  CORRELA  TION  OF  INCREMENTAL  BASE  PRESSURE  DUE  TO  PROPULSI VE 
JET  AS  A  FUNCTION  OF  JET  TO  BASE TOM AX  IMUM  BODY  DIAMETER 
FUNCTION 


Using  correlation  parameters  called  B  and  C  which  involve 
a  reference  base  pressure  level  and  known  information  about 
nozzle  geometry  and  pressure  ratio,  a  remarkable  collapsing 
of  the  available  data  is  obtained,  as  shown  in  Figure  50. 


The  parameters  B  and  C  are  defined  below: 

P  P 

B  =  M  In"1  (.815  -  1.15  InK)  (~)  (~) 

00  P  P 

00  00 

M  K 


where  (g— )  =  model  base  pressure  ratio 

00 

M 


(—•)  =  base  pressure  ratio  at  K  from  Figure  50 


K 


K 


M  Pn  +  M,P. 
®  B  J  D 

- g - 


M.  ®  Mach  number  corresponding  to  nozzle  pressure 
3  ratio 


B 


=  base  perimeter  exposed  to  free-stream 


Pj  =  jet  perimeter 

P  «  perimeter  of  circle  of  area  equal  to  the 
0  base  area 


2  (D./D.)  -  1.5  A 

C  -  ( •  37Mco  >  0.62)  M.  b'-i16 - M.  V 

■^design  -^design  j  AB 


where  D, 


base  diameter 


Dj  =  jet  exit  diameter 

M.  =  nozzle  design  exit  Mach  number 
-'design 


Aj  =  jet  exit  area 
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Figure  50:  ANNULA  R  SUPERSONIC  BASE  PRESSURE  CORRELA  TION 


=  base  area 


ab 

P.  =  nozzle  exit  static  pressure  (just  inside 
-1e  the  nozzle) 

The  parameter  (P^/P^)^  is  plotted  as  a  function  of  k  in 
Figure  51.  B  is  plotted  as  a  function  of  C  in  Figure  50. 
The  procedure  for  determining  base  drag  is  then: 

1 .  Compute  C 

2.  Find  the  corresponding  B  from  Fiqure  49. 

3 .  Compute  K 


4. 


Find 


from  Figure  50. 


A  method  to  account  for  upstream  boundary  layer  effects  has 
been  developed  by  Nash  (Reference  26).  Nash  proposed  a 
curve  of  the  limiting  value  of  the  base  pressure  coefficient 
for  zero  boundary  layer  thickness  as  a  function  of  Mach  number 
Then  he  obtained  an  impressive  correlation  (Figure  52)  for 
the  increment  in  base  pressure  coefficient  due  to  a  finite 
momentum  thickness  0  approaching  a  base  of  effective  height 
h.  Thus  the  base  pressure  coefficient  in  the  absence  of 
jet  effects  can  be  estimated.  The  momentum  thickness,  0,  to 
make  use  of  the  correlation  can  be  obtained  from  experimental 
data,  detailed  theoretical  calculations,  or  engineering 
procedures,  as  time  permits. 


5.2.4  Scrubbing  Drac 


Provisions  for  scrubbing  drag  calculations  follow  the  method 
of  Reference  27,  which  has  been  used  to  correct  the  747  JT9D 
engine  uninstalled  performance  for  the  scrubbing  drag  due  to 
the  fan  nozzle  flow  impinging  on  the  nacelle  afterbody. 
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Figure  52:  CORRELATION  OF  BASE  PRESSURE  DATA  FOR  SUPERSONIC  FLOW 


The  method  divides  the  calculation  of  viscous  scrubbing  drag 
into  two  parts:  The  portion  of  the  flow  from  the  nozzle  exit 
to  the  end  of  the  core  flow,  and  the  portion  from  the  end  of 
the  core  flow  to  the  end  of  the  afterbody.  The  flows  differ 
in  that  the  core  region  has  a  constant  velocity  scrubbing  the 
surface,  and  the  latter  portion  is  in  a  region  of  velocity 
decay  (Figure  53) .  For  the  region  of  the  core  the  viscous 
drag  is  expressed  as 

Dc  =  0.036bpv'2VJ1,8Xc'8 

where  b  =  width 

p  =  density 

v  =  kinematic  viscosity 
Vj  =  exhaust  flow  velocity 

Xc  =  core  length 

For  the  region  beyond  the  core,  the  drag  is  expressed  as 


D  =  0.0288b 

i  \ 

/ 

XC  1 

1  dv  +  c 

J  Y  ^ 

xc 

xc 

„  .25  .7T2T 

X  J 

Where  Vx  is  the  local  velocity  at  station  K. 

The  core  length  (X_)  and  velocity  decay  (V  =  f(x-X  ))at  any 
point  downstream  or  the  core  are  determined  using  the  method 
of  Lawrence  (Reference  27) .  The  maximum  velocity  downstream 
cf  the  core  was  correlated  as 


VMAX 
Vj  -  V 


V 


00 


00 


=  K  (££) 


B 


Where 

VM-v  =  maximum  velocity  downstream  of  X^^ 
MAX  CORE 

Vj  =  fully  expanded  jet  velocity 

Vto  =  free-stream  velocity 


113 


♦  /•  '' 


K  =  velocity  decay  constant 

B  =  velocity  decay  exponent 

AX  =  local  station  downstream  of  X,,-.-.,, 

h  =  flow  height 

The  constants  K  and  B  were  developed  from  wind  tunnel  data, 
and  are  presented  in  Figure  54.  Intermediate  values  for 
other  Mach  numbers  are  assumed  proportional  to  free-stream 
ram  pressure  ratio.  Core  length  is  calculated  by  setting 
Vj  =  and  solving  directly  for  Xc. 

The  above  method  exists  now  in  subroutine  form  (see  Volume  II) . 
5.3  ANGLE-OF -ATTACK  EFFECTS 

Angle-of-attack  effects  on  nozzle  internal  performance  and 
nozzle/afterbody  drag  are  difficult  to  predict  using  techniques 
suitable  for  preliminary  studies  because  the  existing  body  of 
experimental  data  to  show  these  effects  is  very  limited.  In 
addition,  the  data  available  are  for  specific  configurations 
and  cannot  readily  be  generalized.  Reference  31  contains 
experimental  nozzle/afterbody  data  at  angle-of-attack  for  a 
twin-engine  fighter  at  Mach  0.9.  These  data  show  that  there 
is  negligible  effect  of  angles-of-attack  from  -2  to  +  10 
degrees  on  nozzle/afterbody  drag  for  the  tails-off  configuration. 
When  the  tails  are  on,  however,  there  is  a  large  increase  in 
drag  due  to  angle-of-attack.  There  is  not  enough  data,  however, 
to  determine  how  much  the  angle-of-attack  drag  increase  depends 
on  the  particular  configuration  used. 

Variations  in  angle-of-attack  can  affect  the  nozzle  thrust 
coefficient  by  changing  the  local  pressure  field  at  the  nozzle 
exit,  and  by  introducing  asymmetries  and  interactions  into 
the  nozzle  exit  flow  field  due  to  thinning,  thickening,  or 
separating  the  boundary  layer  ahead  of  the  nozzle.  The  extent 
to  which  these  effects  alter  the  nozzle  thrust  depend  on  many 
variables,  including  nozzle  geometry,  pressure  ratio,  afterbody 
shape,  fin  location,  and  angle-of-attack  range.  No  generalized 
data  were  located  to  predict  these  effects.  It  can  be  expected 
however,  that  for  a  choked,  underexpanded  nozzle  at  reasonably 
low  area  ratio,  the  effects  of  moderate  changes  in  angle-of- 
attack  on  nozzle  thrust  coefficient,  Cp  ,  can  be  neglected  for 
purposes  of  preliminary  studies.  g 

For  most  preliminary  design  studies ,  steady  state  cruise 
performance  is  emphasized,  therefore,  the  effect  of  changes 
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in  angle-of-attack  on  afterbody  drag  is  not  of  primary 
interest  and  can  also  be  neglected. 

If  angle-of-attack  performance  is  considered  critical  to 
the  design,  a  more  extensive  study  should  be  made  of  the 
particular  configuration  of  interest  to  determine  whether 
the  effect  is  likely  to  be  large  or  small. 

5.4  INSTALLATION  EFFECTS  OF  THRUST 

VECTOR  CONTROL  AND  THRUST  REVERSERS 

The  installation  of  thrust  vector  control  and  thrust  re  terser 
equipment  (in  the  stowed  position)  can  create  losses  in 
internal  thrust  coefficient,  increases  in  nozzle/afterbody 
drag,  and  increases  in  airplane  weight. 

To  determine  the  magnitude  of  these  losses  and  weight  increases 
it  is  necessary  to  compare  the  performance  and  weight  of 
various  airplane  design  concepts  which  use  thrust/vectoring 
and  reversing  installations  with  similar  airplanes  that  do  not 
have  TR/TV  installations.  Only  a  limited  amount  of  data  are 
available  from  studies  of  this  type,  and  the  configurations 
are  mostly  of  the  V/STOL  type.  The  results  from  recently 
completed  Air  Force  Contract  F33615-71-C-1850,  (Reference  29) 
are  used  to  illustrate  the  magnitude  of  the  effects  involved. 
The  results,  presented  in  Table  IV,  were  obtained  using  a 
subsonic  STOL  transport  as  a  basepoint.  Various  thrust 
vectoring  and  reversing  installations  were  designed  for  the 
transport  and  the  final  performance  and  weights  were  compared 
with  the  basepoint  to  obtain  the  increments  shown  in  Table  IV. 
For  a  full  description  of  the  configuration,  see  Reference  29. 
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TABLE  IV  (Continued) 
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Based  on  the  data  presented  in  Table  IV,  it  appears  that 
the  installation  of  thrust  vectoring  and  reversing  equipment 
can  change  the  nozzle  internal  thrust  coefficient  by  as  much 
as  -0.024  to  +0.024  and  increase  the  weight  from  620  to 
2680  lb.  for  a  typical  subsonic  STOL  airplane.  Since  the 
weight  effect  appears  to  be  the  most  significant  effect,  it 
is  expected  that  the  effect  of  these  installations  in 
supersonic  aircraft  would  result  in  somewhat  more  severe 
penalties  to  the  mission  range  due  to  the  decreased  L/D  ratios 
as  design  Mach  number  is  increased. 
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*T ablet  indicated  by  arabic  numbers  (in  Figures  55  -  59)  such  as 
Table  7,  refer  to  the  tables  sketched  in  Figures  9  and  10. 

PERF  INPUTS:  PT  ,  ,  Me,  WB  ,  Wr  ,  7,  TABLES  1  -  7,  * 

'<»  5  eENG  USEC 

AND  TABLE  8  IF  SCHEDULED-SPILLAGE  MODE  IS  USED 
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1  RETURN  WITH  gp  .  P|NLET 

Figure  55.  FLOW  CHART  FOR  INLET  PERFORMANCE 
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Figure  58:  FLOW  CHART  FOR  SCHEDULED-BYPASS  MODE 
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Figure  61:  FLOW  CHART  FOR  SUBSONIC  AND  TRANSONIC  INLET 
INTERNAL  PERFORMANCE  CALCULATION 
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Figun  62:  FLOW  CHART  FOR  SHOCK  LOSS  CALCULA  TION 


Fifure  65:  FLOW  CHART  FOR  B  YPASS  FLAP  DRAG  CALCULA HON 
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Flight  Dynamics  Laboratory 

Air  Force  Systems  Command 
Wright-Patterson  Air  Force  Base,  OH 

This  report  presents  the  results  of  a  research  program  to  develop  a 
procedure  for  use  in  calculating  propulsion  system  installation  losses. 
These  losses  include  inlet  and  nozzle  internal  losses  and  external  drag 
losses  for  a  wide  variety  of  subsonic  and  supersonic  aircraft  configu¬ 
rations  up  to  Mach  4.5.  The  calculation  procedure,  which  was  largely 
developed  from  existing  engineering  procedures  and  experimental  data,  is 
suitable  for  preliminary  studies  of  advanced  aircraft  configurations. 
Engineering  descriptions,  equations,  and  flow  charts  are  provided  to 
help  in  adapting  the  calculation  procedures  to  digital  computer  routines! 
Many  of  the  calculation  procedures  have  already  been  programmed  on  the 
CDC  6600  computer.  Program  listings  and  flow  charts  are  provided  for 
the  calculation  procedures  that  have  been  programmed.  The  work  accom¬ 
plished  during  the  program  is  contained  in  four  separate  volumes. 

Volume  I  contains  an  engineering  description  of  the  calculation 
procedures.  Volume  II  is  a  programmers  manual  containing  flow  charts, 
listings,  and  subroutine  descriptions.  Volume  III  contains  sample 
calculations  and  sample  input  data.  Volume  IV  contains  bookkeeping 
definitions  and  data  correlations. 
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